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MANY-BODY THEORY OF A RAPIDLY VARYING, 
INHOMOGENEOUS ELECTRON GAS 


By John William. Gadzuk 
Electronics Research Center 

SUMMARY 

The case of an inhomogeneous electron gas within which the density 
variation is significant over a spatial range of the order of a Fermi wave- 
length is considered in this report. It is seen that for most systems of 
physical interest this sort of non- uniformity is a result of diffraction effects. 
This is a fundamentally different phenomenon than can reasonably be treated 
by the density gradient method of Kohn for slowly varying inhomogeneous 
electron gases. Several sample cases are treated. 

The first considerations are directed towards the problem of a weak 
periodic potential in an interacting electron gas. The momentum -dependent 
self- energy is calculated for an electron propagating in the many-body me- 
dium of an electron gas plus a periodic lattice pseudo-potential. This is the 
equivalent of a quasi-particle energy spectrum and thus an orthogonalized 
plane wave energy band. It does not appear that the lattice drastically 
changes qualitative aspects of plane wave many-body theory. 

A dielectric formulation for a general inhomogeneous electron gas is 
presented. By introducing a new image technique, the dielectric function 
within the random phase approximation, which is valid in the surface region 
of an electron gas, is obtained. 

A Green's function formalism is developed for treating the static dielec- 
tric screening of a point impurity in an electron gas. 

The surface dielectric function is used with the impurity screening 
formalism to treat the problem of impurity screening in the surface region. 
This is an idealized model of ionic adsorption on metal surfaces. Screening 
charge densities resulting from volume polarization effects are calculated. 
From these results, it is seen why unjustifiable application of classical 
image forces in previous adsorption theories has fortunately produced rea- 
sonable results. 

A new method for obtaining the appropriate plasmon contribution to the 
electron self- energy in the surface region is developed. With these results, 
the electron gas surface potentials calculated by Loucks and Cutler are then 
improved. 



I. INTRODUCTION 


The striking success of the single - electron theory of metals has pro- 
vided large amounts of quantitative information pertaining to electron 
behavior in real metals. The initial philosophy upon which the single- 
electron theory was based concerns an electron in an N-electron metal 
feeling a static periodic potential provided by the ionic cores of the nuclei 
and inner-shell electrons, plus the time-averaged, self-consistent poten- 
tial of all other N-l conduction electrons. The resulting energy eigenfunc- 
tions for such a system are Bloch functions which reflect the symmetry of 
the particular lattice in question. The fact that the pseudo-momentum of an 
electron is a good quantum number solely as a result of the particular trans- 
lational invariance present is made manifest in the Bloch function represen- 
tation. In a periodic potential, the translational invariance is that of the 
translation group of the crystal only. Thus in the reciprocal lattice of the 
crystal, momentum is a good quantum number only within an additive recip- 
rocal lattice vector and, consequently, crystal momentum is the conserved 
quantity usually as sociated with translational invariance. This situation 
has given rise to an energy band theory in the reduced zone scheme and has 
required the introduction of a variable called the band index to keep track of 
the fact that crystal momentum and not mechanical momentum is conserved. 
The Fermion nature of an electron is considered in this picture by allowing 
only two electrons to occupy any Bloch state, the factor of 2 arising from the 
two different eigenstates of spin. The ground state of the system is obtained 
by having electrons occupy the N/2 lowest Bloch states. In this picture, 
which is called the Hartree approximation, there are absolutely no correla- 
tions between electrons either of a dynamically or statistically imposed 
nature. The original N-particle Hamiltonian has been decoupled and reduced 
to a sum of N- single-particle Hamiltonians. The wave function for the sys- 
tem is a product of N-non-interacting wave functions. 

The next advance in single- electron theory was the inclusion of exchange 
effects. It was realized that some correlation existed between electrons of 
parallel spin as a result of the exclusion principle. Electrons in the same 
spin state tend to stay away from each other. Thus, when evaluating the 
time-average, single- electron interaction with all other N-l electrons, 
there is a deficiency of parallel spin electrons around the given electron. 

This means that the repulsive interaction is lessened and thus the total 
energy of the system is reduced, in fact, by a sufficient amount to allow for 
binding, although not by the quantitatively correct amount. Decoupling the 
N-particle Hamiltonian in this way and writing the full wave-function as a 
Slater determinant is known as the Hartree-Fock approximation. 

These single-particle approximations have been quite useful in gaining 
much understanding of the behavior of metals. Single- electron energy band 
theory has for the most part successfully explained the characteristic optical 
spectrum of solids. Fermi surface experiments are in relatively good 
agreement with single- electron theory. Single-electron transport theory in 
solids has been successful. And so the list could go on. 
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However, there is a wide class of phenomena which is either not in 
agreement with single- electron theory or categorically has no hope of being 
interpreted in terms of single-electron theory. Examples of disparity 
between single- electron theory and experiment can be found in some cal- 
culated density of states at the Fermi surface, electron specific heats, 
electron effective masses, and bandwidths. Phenomena which cannot be 
treated within the context of single- electron theory include total energies 
of metals, all plasma effects, screening effects, and many electron surface 
effects. It is this class of phenomena which are treated in this document, 
the many-body effects, in hopes that an effective coupling between single- 
electron theory and many-body theory can provide adequate representation 
of nature. 

Wigner (ref. 1) was the first person to offer a semi- quantitative analysis 
of dynamical correlations between electrons of anti-parallel spin. He 
realized that electron- electron repulsion tends to keep anti-parallel spin 
electrons apart so when one calculates total energies of a metallic or elec- 
tron gas system, one should use wavefunctions which consider this correla- 
tion effect. Since the exact wavefunction keeps all electrons apart, they do 
not feel their mutual repulsion so strongly as when in a configuration 
corresponding to uncorrelated wavefunctions. This reduction in energy, the 
difference between the exact ground- state energy and the Hartree-Fock 
ground- state energy, he called the correlation energy. Within an approxi- 
mation which becomes exact in an unphysical, low-electron-density region, 
Wigner then made probably the first many-body calculation. 

The problem of electron correlation and many-body effects then lay 
dormant for many years until the landmark work of Bohm and Pines (ref. 2). 
They modelled a metal in the following manner. Assume the periodic poten- 
tial to be sufficiently weak, such that an adequate first approximation is to 
assume the ion core charge density uniformly smeared over the unit cell. 

The resulting picture is that of an electron gas sloshing around in a uniform 
potential in a manner quite analogous to a classical plasma. This has often 
been referred to as jellium. To consider e-e interactions proceed as 
follows. The e-e Coulomb repulsion is written as a Fourier series. The 
q - 0 (infinite wavelength) Fourier component of the repulsion exactly 
cancels the infinite wavelength attraction of the smeared- out ion charge. 
Problems arise because the unscreened Coulomb l/r interaction is long- 
range. Thus the low q components of the e-e interactions give rise to 
divergences in momentum integrals. Bohm and Pines realized that, in 
truth, the e-e interaction is really screened and, beyond some character- 
istic length of the order of a Thomas-Fermi screening length, an electron 
reacts with other electrons in some average way. They envisioned the 
electron gas sloshing around as a cooperative phenomenon which could be 
expressed in terms of a Fourier sum of density fluctuations about the mean 
electronic density. These density fluctuations were just the plasma oscilla- 
tions of the electron gas. Since the plasma oscillations extend over a large 
region of space, the low-q components of the e-e interaction could be seen 
to affect the long wavelength density fluctuations. Thus, if one could sep- 
arate the full interacting electron Hamiltonian into a plasma part and an 
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electron part, then the divergent-appearing, low-q portion of the e-e inter- 
action would go into theplasma Hamiltonian and no divergence problems 
would appear in the quasi-electron part which uses the high-q portion of the 
unscreened interaction. To carry out this procedure, they performed a 
series of appropriate canonical transformations on the interacting electron 
Hamiltonian and obtained a Hamiltonian which was conveniently separable 
into a plasma part, an electron-like part, and a weak coupling between the 
two which could usually be treated in perturbation theory. The plasma part 
was just a sum of harmonic oscillators which could then be quantized in the 
appropriate manner for a Bose field and the resulting quanta were then 
called plasmons. The characteristic plasma frequency for electrons of a 
density found in metal conduction bands was seen to be of the order of 
10 eV, much higher in energy than could be excited by e-plasmon interactions 
or thermal effects. Thus the only effect of the low-q portion of the e-e 
interaction on the total energy of the system was the addition of a plasmon 
zero point energy l/2 h top. In any event, the recognition of the electron gas 
as a field in which density fluctuations could occur and the subsequent recog- 
nition of the many-body cooperative phenomena giving rise to plasmons 
spurred a wave of new work on electron correlation effects in electron gases 
and real metals. It should be noted that, in the Bohm-Pines (BP) theory, 
first use of the approximation, which has been the blessing and nemesis of 
many-body theorists, the random phase approximation (RPA), was made. 
Unfortunately this approximation, the nature of which will be discussed later, 
makes the BP and most many-body theories rigorously valid only for elec- 
tron gas densities much higher than those found in conduction bands. 

The pioneering work of Bohm and Pines (BP) led Gell-Mann and 
Brueckner (G-MB) to the most significant breakthrough in modern many-body 
theory (ref. 3). The idea of the electron gas as a field in which quantized 
density fluctuations could be induced by a particle propagating in this field was 
totally reminiscent of the non- relativistic specialization of the field theoretic 
ideas developed in quantum electrodynamics. In quantum electrodynamics, 
success in the usage of perturbation series generated in a propagator form- 
alism and pictorially represented with Feynman diagrams had been realized. 
What was the difference between an electron, which in physical reality can- 
not be imagined as decoupled from the electromagnetic field, propagating in 
the electromagnetic field and inducing density fluctuations and an electron 
coupled to the field of the many-body medium within which it propagates and 
induces fluctuations? Mathematically, as well as physically, the distinction 
was small. For the Fermi gas, the following picture thus evolved. Initially, 
with no e-e interactions, we can define a Fermi sphere in momentum space 
which is the equivalent of the field theoretic vacuum state. Turning on e-e 
interactions results in the excitation of two high-momentum propagating 
electrons above the Fermi sphere and two low-momentum holes propagating 
beneath the Fermi sphere. The electron-hole pairs can then annihilate, 
create new e-h pairs, and undergo complicated exchange processes, or a 
whole host of various processes. Gell-Mann and Brueckner set up an 
infinite order perturbation series to describe the effects of the e-e inter- 
action. Each term in the perturbation expansion could be represented by 
a Feynman diagram. They reckoned that one could approximate the exact 
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series by considering only the most divergent diagram in any order of the 
perturbation expansion. Through this assumption they were able to choose 
a particular set of diagrams, the bubble or ring diagrams to be discussed 
later, the infinite series of which was exactly summable. Thus the stigma 
of a perturbation theory treatment was removed. As it turned out, by 
choosing the particular, most divergent diagrams that they did, they ob- 
tained exactly the BP result for the correlation energy. Hence, it can be 
seen that the RPA is made by choosing the ring diagrams only. It is within 
the context of a diagrammatic analysis that we chose to understand and use 
the RPA. One limiting feature of the G-MB theory was that it did not ex- 
plicitly display the existence of the plasma poles. However, the theory did 
open the road for the vast usage of field theoretic and Green’s function tech- 
niques which have proved to be the forte of subsequent many-body theory. 

Sawada (ref. 4) and Sawada, Brueckner, Fukuda, and Brout (ref. 5) pro- 
vided an elegant formulation of the G-MB theory in terms of the equation- of - 
motion method. This work was stimulated by Sawada et al who realized the 
mathematical similarity between the G-MB electron gas pair theory and the 
theory of scaler-pair meson fields. In their formulation of the electron gas 
problem they showed that the G-MB theory, in fact, had the plasma pole 
implicit within it. 

DuBois (ref. 6) presented extensive work on the formal field theoretic 
aspects of the electron gas. He showed that great and almost impossible 
effort was required to extend most aspects of the RPA to the region of 
metallic density when done through a systematic application of perturbation 
theory. Much work was done on the plasmon aspects of the electron gas. 

The connection between the propagator for a density fluctuation and the wave 
number and frequency-dependent dielectric function was given. Thus, if 
one can calculate a dielectric function or equivalently a polarizability for 
the electron gas, then the formal apparatus for calculating other physically 
relevant quantities, such as self-energies, correlation energies, and plas- 
mon damping "constants," exists. 

Another tack was taken by Hubbard (ref. 7) as he tried, in an approxi- 
mate manner, to make the RPA applicable for electron gases of metallic 
density. For many properties, the RPA and Hubbard-corrected RPA give 
quite similar numerical results. In any event, Hubbard’s work seems to be 
the best correction to the RPA presently available. 

Ehrenreich and Cohen (ref. 8) presented an elegant piece of work which 
demonstrated the equivalence between a linearized, self-consistent field 
theory and the RPA with Hartree factorization of two particle operators. 

They obtain simple expressions for the dielectric functions. 

More recently Hedin (ref. 9) considered the electron gas from a slightly 
different point of view. A Green's function theory is developed in which per- 
turbation expansions are done with a screened rather than bare Coulomb 
interaction. As expected, this greatly improves convergence problems. He 
also provides large quantities of numerical tabulations of relevant electron 
gas properties. 
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Hedin, Lundqvist, and Lundqvist (ref. 10) have used a novel approach 
to evaluate electron and hole self- energies resulting from a coupling between 
virtual electrons or holes and virtual plasmons. Although DuBois (ref. 6) 
warned that this procedure has no physical validity, the author of reference 
10 seemed to show the existence of a new, many-body resonance, a coupled 
hole-plasmon state, which may have been observed and which, if it in fact 
exists, explains a sizable quantity of anomalous density of state data for 
such materials as nickel. 

The preceding review of published work has been given to set the stage 
for the problems of central interest in the present thesis. In all of the 
aforementioned theories, the jellium model has been treated; that is, an 
electron gas in which complete translational invariance exists. This 
approach is the reasonable first attack on the problem of many-body effects 
in metals. However, in real metals there are many different occurrences 
which destroy the assumed homogeneity of space. The obvious first observa- 
tion is that the ion core potential is not smeared out over the unit cell; a 
periodic potential is present instead. This aspect of nature should be con- 
sidered before the final many-body theory answer is realized. A second 
class of phenomena concerns the existence of impurities or imperfections in 
the metal or electron gas. The third major class of problems might be re- 
garded as a subclass of the imperfection group, problems arising from the 
presence of a surface or boundary in an electron gas. It is these three types 
of problems which will be treated here. As will be seen in the historical re- 
view and subsequent chapters, there are many mathematical and physical 
problems which arise in this type of many-body problem and, for the most 
part, both the theory and mathematical apparatus for dealing with these 
problems remain in their infancy. 

Some brief considerations of the work which has been done with the 
lattice in an electron gas are now presented. 


In the self-consistent field theory of Ehrenreich and Cohen (E-C) (ref. 8), 
they consider lattice effects in an electron gas by doing expansions in a Bloch 
wave rather than plane wave basis. They obtain a dielectric function which 
can be split up into terms involving intra-band and inter-band transitions 
from which suitable approximations on each part can be made. 

Falk (ref. 11) has presented a calculation on the influence of a lattice 
upon the dielectric properties of an electron gas. Again, expansions are 
done in a Bloch-wave basis. By considering the most divergent portions of 
integrals, he obtains relatively simple formal expressions for dielectric 
functions and correlation energies. 

Further work along these lines was done by Pratt (ref. 12) in an attempt 
to introduce self-energy effects into band theory. He obtains single-particle 
Schrodinger-type equations which include periodicity and many-body effects 
but which are almost impossible to solve self- consistently , at least from the 
practical point of view. 
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The most well known attempt to put a lattice into an electron gas is the 
projected wave-field approach of Bassani, Robinson, Goodman, and 
Schrieffer (ref. 13). Essentially what they do is the following. The exact 
wave functions are expanded in a few orthogonalized plane waves (OPW). 

They then define two electron fields, one corresponding to the plane wave por- 
tion of the OPW's and the other corresponding to the projection onto the core 
states. Each field is quantized independent of the other and a weak coupling 
between the two fields might then exist. From this point of view they obtain 
formal energy spectra for quasi-particles which contain both periodicity and 
self- energy effects. 

Sham and Kohn (ref. 14) made reference to lattice effects when they 
presented one of their works on inhomogeneous electron gases. Although 
their theory is quite formal, applicable mathematical techniques are pro- 
vided. 

The author (ref. 15) has investigated this problem from a point of view 
which in many ways is similar to that of the authors of reference 13. How- 
ever, the mathematical aspects are carried further towards completion and 
contact with physical pictures is maintained. This is discussed in detail in 
section III of this document. 

The second interesting class of translational invariance spoilers are 
impurities in an otherwise homogeneous electron gas. The simplest example 
of an impurity in such a gas is an atom. If one assumes a liberal point of 
view, it can be seen that an atom is nothing but an electron gas with a point 
impurity immersed in it -- the nucleus. Although nuclear physicists re- 
ject the idea that their life's work is concerned with the study of a piece of 
"dirt, " for present purposes this is a sufficiently acceptable outlook. The 
electron gas tries to screen the nucleus, Friedel oscillations are set up, 
and thus shell structures result, or so it was thought. 

This was the point of view taken by Baraff and Borowitz (ref. 16) and 
DuBois and Kivelson (ref. 17) as they tried to improve upon Thomas-Fermi 
atoms by systematically using many-body theory ideas and density gradient 
expansions in hopes of generating shell structure. As they learned, for 
subtle reasons this is a theoretically impossible task. 

In a totally new approach to the inhomogeneous electron gas problem, 
Kohn and Sham (ref. 18) were able to obtain almost perfect shell structure 
radial charge densities for atoms when treated as a many-body screening 
problem. 

The next class of impurity problems is that of a point impurity 
immersed in a metal. For the most part, the case of substitutional rather 
than interstitial impurities is considered in this document. Langer and 
Yosko (ref. 19 ) presented the first calculation on many-body impurity 
screening effects. Essentially they considered the impurity to be a source 
which interacts with the electron field, thus causing an energy shift in single- 
particle energies. This energy shift is then related to screening charge 
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densities in a self-consistent manner through Poisson's equation. They pro- 
vide extensive numerical computations of screening charge densities calcul- 
ated within the RPA and Hubbard approximation. As expected in a correct 
treatment, long-range Friedel quantum density oscillations are found which 
decay in amplitude as r“^ and oscillate with a spatial frequency of twice the 
Fermi momentum. This result will be important in later sections of this 
thesis. 


A purely Green's function treatment of point Coulomb impurities was 
offered by Eayzer (ref. 20). He shows the relationship between the purely 
diagrammatic self-consistent field approach of Langer and Vosko and the 
Green's function method. 


An extremely novel method of tackling the impurity problem came from 
Sziklas (ref. 21). He combined the Bohm-Pines plasma theory with the E-C 
self-consistent field theory to show that, under certain conditions, the 
screening charge around an impurity could be regarded as a bound state 
plasmon, a trapped density fluctuation. The primarily s- state plasmon 
would be expected to have a discrete hydrogenic-like excitation spectrum 
which in principle would be observable. It is a shame that this idea has not 
been further developed. 

The more traditional approach to the impurity problem in a metal is 
based in large upon the work of Friedel (ref. 22). If screening is regarded 
as a scattering process, describable in a phase shift analysis, then his 
famous sum rule gives the required relation between various phase shifts 
to guarantee self-consistency in the screening charge density. 

Although the tight- binding impurity theory of Koster and Slater (ref. 23) 
does not on the surface appear to use many-body screening effects as a 
basic ingredient, in fact it does so implicitly. By considering an impurity 

so localized that it only gives non- vanishing matrix elements between 
Wannier functions centered on the impurity site, it is implicitly required 
that the impurity be rapidly screened. This formalism has enabled Sokol off 
(ref. 24) to treat such problems as magnetic impurities in terms of band 
wave -functions following the lines of Anderson's theory (ref. 25). By using 
the various matrix elements as adjustable parameters, many-body screen- 
ing effects are introduced phenomenologically. 

Clogston (ref. 26) presented a method of relating the Koster-Slater 
matrix elements to the Friedel sum rule such that self-consistency was ob- 
tained. 

Recently Sham (ref. 27) has attacked the impurity problem from a many- 
body theory point of view with a direct application to semi-conductor or in- 
sulator impurities. However, most of the formalism is quite useful in the 
metal impurity. 

The author has considered the metal impurity problem from the point of 
view outlined in the Thomas-Fermi atom theories (ref. 28). He discovered 
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several hidden assumptions which are required but rarely stated in order to 
obtain the standard RPA screening charge densities. Certain aspects of this 
work are discussed in section IV of this technical report. 

By far, the most significant work on the problem of inhomogeneous 
electron gases has emananted from Kohn and others working with him, 
namely, Hohenberg, Sham, and Majumdar (refs. 29-31). The original basis 
of this work is found in a presumably general theory which states that all 
properties of an inhomogeneous electron gas can be expressed in terms of a 
unique, universal functional of the electron density at a point (ref. 29). Of 
course no method has yet been devised to obtain this functional, but the 
theory has given cause for more work using density as the basic variable. 
Thus slowly varying or small departures from the mean electron density 
are treated in a formalism based upon perturbation-like expansions in 
powers and orders of the density gradient. In general, polarizabilities or 
dielectric functions calculated for the uniform gas are used in some sort of 
perturbation treatment. Another innovation has been to use a WKB basis as 
a starting point for expressing non-uniform potentials and densities (ref. 30). 
Good preliminary success has been obtained. 

Kohn and Majumdar also investigated the conditions required for an im- 
purity to be screened by formation of a localized bound state pulled out of 
the band of continuum states or by formation of a virtual bound state within 
the occupied portion of the band. Although striking in some of its successes 
in treating slowly varying inhomogeneities, it is the author's feeling that 
this approach cannot be justifiably extended to the cases of lattice or surface 
effects for the following reasons, which are similar to the reasons for the 
failure of the many-body Thomas-Fermi atom theory. There are two 
fundamentally different origins of density variations in an electron gas. In 
the case of a slowly varying potential, the number of non-interacting 
electrons with positive kinetic energies roughly follows the contours of the 
potential variation. In this case it is possible to make a one-to-one 
correspondence between the total energy and the density. This is the type 
of situation for which the HK theory is intended (ref. 29). On the other hand, 
for the lattice- or surface-induced inhomogeneities, the density variations 
are caused by interferences set up between electrons when each non-inter- 
acting electron is forced to satisfy the Bloch or surface boundary condition, 
or, in other words, phase correlations give rise to rapidly varying total 
density changes. In this case, there is no exclusion of low kinetic energy 
electrons in the region of low total density so a one-to-one correspondence 
between density and energy does not necessarily exist. If this conjecture is 
valid, then the "density as a basic variable theory" cannot, in its present 
form, be expected to be applicable to lattice- or surface-induced inhomo- 
geneities. It should be pointed out in passing that this is precisely why a 
systematic usage of gradient expansion techniques failed in the atom theories, 
for in atoms the Friedel oscillations or shell structure occur basically as a 
result of forcing a boundary condition on the non-interacting electron wave- 
functions, that of having zero amplitude at the origin. This is not a real 
many-body effect but instead a property of the wave nature of electrons which 
would be impossible to describe in terms of many-body interactions. Further 
considerations along these lines should be forthcoming. 
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The methods of Kohn have been further carried on by Grover (ref. 32). 
Essentially he does a density matrix formulation using a 'WKB basis to sym- 
bolize the inhomogeneous nature of the system. 

The third class of inhomogeneity effects in an electron gas, and probably 
the most neglected in spite of their importance, is that of surface effects. 
Perhaps the plethora of phenomenological surface theories have led people to 
conclude erroneously that surface theory is progressing at a satisfactory 
rate. One must be careful not to fall into the tempting trap of believing that 
reproduction of numbers in agreement with experiment constitutes under- 
standing. With a few exceptions, this has been the curse of surface physics 
and has led to its retarded theoretical growth. 

An amazing feature pertaining to surface effects is that, almost in toto, 
electronic surface properties are manifestations of pure many-body effects. 
For instance, it is well known that the surface dipole barrier contributes 
only a very small percentage, say 10 percent, to the surface barrier, that 
barrier which allows for the binding of a metal. If we assume an idealized 
model with no surface dipole, then in the Hartree approximation, where any 
attractive potential would be electrostatic, a metal does not exist. Conse- 
quently, the total surface barrier arises from exchange and correlation 
effects which are of purely quantum mechanical origin. On the metal side 
of the surface, an electron is propagating through the many-body field within 
which it can undergo exchange processes or induce density fluctuations and 
thus self-energy processes resulting in a lowering of its energy whereas , on 
the vacuum side, the electron is completely decoupled from such a field and 
thus is found to be at a higher energy. Any "theory" which does not recog- 
nize this picture of a surface is categorically in conflict with standard 
principles of physics as stated through Maxwell's equations, the Pauli Ex- 
clusion principle, the correspondence principle, and the principle of com- 
plementarity. 

With regards to the surface problem, one of the immediately relevant 
objects of interest is knowledge as to the form of the average total potential 
(exchange and correlation) as one approaches the surface from deep within 
the metal. The first work along these lines was done by Bardeen (ref. 33). 

He showed that using a sine wave basis with a mathematical nodal surface 
somewhere outside the physical surface was a good approximation to the 
traditionally desired basis consisting of decaying exponentials outside the 
physical surface. Es sentially the justification is as follows. The sine 
waves depart from decaying exponential behavior only in the region where the 
wave function amplitude becomes very small. Since it really doesn't matter 
how nothing is described, the most convenient way available should be used 
as resulting errors will be negligible. With this in mind, Bardeen calculated 
the values for the exchange energy at two points in the surface region, and 
then formed a smoothly oscillating curve which merged into the image 
potential far outside the metal and which joined onto the interior potential 
deep within the metal. 
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Several years later, Juretschke (ref. 34) performed numerical calcula- 
tions of the exchange integral for a continuous distribution of points in the 
surface region. The resulting potential approached the interior potential 
within the metal and oscillated with increasing amplitude as the surface was 
approached, thus giving a functional form for the surface barrier. An 
especially clear physical picture of the origin of the reduced exchange energy 
attraction in the surface region in terms of squashed exchange holes is 
given. 


Loucks and Cutler (ref. 35) added further sophistications to the 
Juretschke picture by realizing that a screened Coulomb exchange potential 
should be used. Furthermore they included correlation effects in an approxi- 
mate manner within the context of the BP plasma oscillation theory. Further 
work along this line has been completed by the author and constitutes section 
VII of this document. 

One of the more dramatic and consequently interesting electronic sur- 
face phenomena is the existence of a collective mode of oscillation with a 
characteristic frequency of the bulk plasma frequency divided by the square 
root of 2. Most of the work which has been done on surface plasmons stems 
from the original work of Ritchie (ref. 36) and Stern and Ferrell (ref. 37). 
Very simple derivations and explanations of the surface plasmon are given. 
Essentially the surface plasmon results from the collective oscillations of 
the induced polarization charge in the surface region which is screening out 
the fact that there is a surface. Since the effective electron density in the 
surface region is less than in the interior, it is to be expected that the rele- 
vant plasma frequency would be reduced. 


Ritchie and Marusak have presented detailed calculations of the surface 
plasmon dispersion relation (ref. 38). 

In section IV of this work, it is shown that in a calculation of the dielec- 
tric function appropriate to the surface region of an electron gas, there 
exists a pole at precisely the surface plasma frequency. Of course, this is 
required in a proper treatment. 

The existence of the surface plasmon should be of more than peripheral 
interest to solid-state theorists. There have been several cases in which 
optical properties of solids were obtained by a Kramers-Kronig analysis of 
optical reflectance data which disagreed with photo- emission data. Assign- 
ments of certain interband transitions have been given as a result of the 
reflectance data, only to find out later that the strong absorption was caused 
by the excitation of a surface plasmon and not an interband transition. Such 
transitions, of course, were not observed in the photo- emis sion studies. 

Fedders (ref. 39) has produced a paper in which he considers electron 
gas surface effects from a dielectric and Green's function point of view. 

This work seems to be a big step in the right direction, although it occas- 
sionally dips into the realm of notational obscurity. 



The last type of surface effect to be considered is dielectric screening 
both of a dynamic and static nature. Heine (ref. 40) has suggested the im- 
portance of the role of screening in Auger- effect experiments at a solid 
surface. 

The author (ref. 41) has done a simplified theory of hydrogen-like atom 
adsorption on electron gas surfaces to illustrate the importance of static 
screening of surface impurities on observable properties of the composite 
system. 

The case of the polarizability of a two-dimensional electron gas, such 
as is found in the idealized model of an n-type inversion layer of a Si-Si O^- 
metal system, has been investigated by Stern (ref. 42). He then calculates 
the two-dimensional screened potential around an impurity and estimates 
charge densities which are not in conflict with those obtained by the author 
(ref. 43). Further elucidation of the author's work in this area forms the 
basis of section VI. 

Finally mention should be made of the work of Bennett and Duke (ref. 44) 
on bimetallic interface phenomena. Based on much of the inhomogeneous 
electron gas work cited in this introduction, they have attempted to obtain a 
description of the resulting electronic configuration, from the dielectric 
point of view, when two different dielectric materials are placed in intimate 
contact with one another. 

With the recording of this necessary background material, the intent of 
this report can now be legitimately stated. Basically the general effects in 
an electron gas with limited translational invariance constitute the subject of 
interest. It is felt that the best way to handle this problem is through 
several specific examples rather than a general formalism. That this is the 
point of view adopted is motivated by the feeling that each specific problem, 
be it the lattice, impurity, or surface problem, has its own unique individ- 
ualistic traits which allow for different simplifying mathematical assump- 
tions, while still remaining within the framework of a single field theory. 
However, in no way is it meant to imply that each problem is distinctly dif- 
ferent from the physical point of view. This, in fact, is a major point, viz. , 
that physics is physics and the artificial barriers set up to distinguish band 
theory from many-body theory from impurity theory from surface physics 
are, in fact, artificial. Nature does not distinguish between fields so why 
should one who is trying to describe nature make such a separation? The 
report is best summarized by making recourse to a descriptive outline of the 
structure of the work. 

Since the general mathematical structure of the entire work is to be done 
within the context of a field theoretic or Green's function point of view, per- 
haps the most general point of view presently available, section II, is 
devoted to generalities. A summary of how the quantized field notion follows 
from the interacting Schrodinger equation is given. The concept of Green's 
functions as particle propagators is developed and the relevance of the 
Green's function to physical observables is noted. The notation used through- 
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out the report is established. The self-energy of a particle propagating in a 
polarizable field is then discussed. 

The first problem to be treated is the effect of a weak periodic potential 
in an electron gas. In section III, the self- energy of an electron in the many- 
body medium of the electron gas, plus periodic potential, is obtained. The 
resulting energy spectrum is the equivalent of an OPW energy band, self- 
consistently corrected for many-body effects. 

The central importance of the wave number and frequency-dependent 
dielectric function to many-body effects is recognized in section IV. Follow- 
ing the self-consistent field theory derivation, a dielectric function for a 
general, non-uniform electron gas whose non-interacting single electron 
wave functions can be expanded in an arbitrary sum of plane waves is derived. 
In keeping with the spirit of the self-consistent field theory, the linearized 
version of this expression is retained. A new image technique which enables 
calculation of a dielectric function in the surface region of a bounded electron 
gas from the derived general expression is developed. It is interesting to 
note that the inverse dielectric function has a pole at the surface plasmon 
frequency as it should. 

The theory of impurity screening in the interior of the electron gas, 
based upon the Thomas -Fermi atom theories, is developed in section V. It 
is shown that to obtain the standard "correct" results, several hidden 
assumptions must be made. 


Section VI is concerned with the static dielectric screening of a point 
impurity at the surface. This is of interest for this is precisely that which 
occurs when ions are adsorbed on metal surfaces. Screening charge den- 
sities are calculated, the results of which show why some of the phenom- 
enological theories of adsorption have achieved as good successes as they 
have. Suggestions for new surface experiments are then offered. 

A new approach to the many-body aspects of the surface potential is then 
presented in section VII. The new potentials are compared with previous 
calculations and comments are made. 

Finally in section VIII, general conclusions of this study are made. 


13 



II. FIELD THEORY AND GREEN'S FUNCTION FORMALISM 


At the risk of belaboring the formal aspects of this report, a brief re- 
view of the Green's function formalism and a statement of the notation used 
throughout the document are included at this point. 

For the interested reader, the book by Schultz is by far one of the best 
of the currently available introductory texts to many-body theory (ref. 45). 
Nozieres's book (ref. 46) also provides a good introduction, but really 
serves its most important function to the "professional" many-body theorist. 
Also worthy of mention is the excellent book by Kadanoff and Baym which 
gives a "layman's" version of Martin- Schwinger theory plus much more 
(ref. 47). 

It might at first be asked: Why do things within a field theoretic context? 

The answer is this: The central feature of the types of processes being con- 

sidered is the method in which a particle interacts with a density fluctuation 
in the many-body medium which is induced as a result of the presence of the 
particle. Since this is a natural field theoretic problem, the motivation is 
clear. 


The usual prescription for obtaining a quantized field theory from a 
particle theory is as follows. The particle equation, in this case, the 
Schrodinger equation, is taken to be the field equation. Through a varia- 
tional procedure one obtains a Lagrangian density. The canonically con- 
jugate momenta to the field are given by relations of the form: 


it (x,t) 


8 L {jK h 
8 (*, t) 


Thus, a field Hamiltonian through a Legendre transformation of the 
Lagrangian can be defined. Quantum mechanics is introduced by requiring 
that the field variable and momenta satisfy commutation or anti- commuta- 
tion relations appropriate to the spin of the field. It might be expected that 
if the field equation were written as an inhomogeneous differential equation, 
then a Green's function method would be used to provide formal solutions. 


How the equations of motion for the electron field Green's functions 
follow from a statement of the Schrodinger equation for a single electron 
coupled to all other electrons through two body interactions becomes a matter 
of interest. The time -dependent Schrodinger equation is written: 


( i It + fv*' =/ di 


r 2 v 


r l" r 2 


) ^ + (r' 2 ,t)^ (r^,t)tf> (r^t) (1) 


In Eq. (1), v (rj -rg) is the generalized interaction term which includes e-e 
interactions and could include single- electron interactions with a source. 
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such, as in the impurity problem, if one so desires. $ > is the adjoint of $ . 
Usually the system of natural units in which h. = c = m = 1 is used, the stand- 
ard unit system being restored only when necessary for computational pur- 
poses. Unit quantization volume will also be used. Occasionally the four- 
dimensional notation (r^;t^) = (1) in configuration space and (pjco) = (p) = (p~, p ) 
in momentum space is used. 

To transcribe Eq. (1) into a quantized field equation for a spin 1/2 field, 
it is required that the Ip ’s, now to be regarded as field creation and annihila- 
tion operators, satisfy the anti- commutation relations: 

( 1 ), ( 1 ’)} = 6 (4) ( 1 - 1 ') ( 2 ) 

(1), ^ d’) } =0 

+ (1), ip + (I')} = °’ 

in which ^ (1) destroys a particle at the space-time point (r^,t^) and $^(1) 
creates a particle at fr^,t^). 


Various Green's functions are defined in terms of these single-particle 
operators and the exact, interacting ground state of the many-body system. 
For example, the single particle Green's function is defined as 


Gj ( r"i tfr' , t* ) 



(3) 


and has the following meaning. It is a measure of the probability that if an 
electron is created in the interacting system at (r';F) it will later be found 
at (rjt). The operator T is called a time-ordering operator which requires 
the earlier of the two times t, t' to be, at the extreme right. Essentially the 
expectation value of the operator ipip ' is taken with respect to the exact, 
time -dependent ground state. In a similar manner, two-particle Green's 
functions are defined as 


G 2 ( r l t l» r 2 t 2’ T 1 t i' r 2 t 2^ 


T^^lV ^ (r 2 t 2 ) ^ ^ t (r' 1 t' 1 ) 


(4) 


Now it can be seen how to obtain a differential equation for a single-particle 
Green’s function from Eq. (1). If Eq. (1) is multiplied on the right by ip t 
(r'j, t' ) and the expectation value of the time-ordered product is taken with 
respect to the exact ground state of the system, then the result is: 
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( T [< 1 $i + t- 2 ), ''<v ‘ i>]> 

= (t[J d 3 r z vlTj -■? ’ 2 ) lp + (r" 2 . t 2 + )if(r‘ 2 , t 2 )if frj, ^ _ 

which by virtue of the fact that T commutes with the purely spacelike opera- 
tors, can be written as 



An ambiguity arises with respect to the first expression because the time 
derivative also plays a role in the time ordering, giving rise to a singular 
point when 1=1'. If the time ordering is carried out using the anti- 
commutation of Eq. (2), it can be seen that 

* (vV^Vi, tp])= ( 6 ) 

- 6 (3) (r-j -rj) 6 (t 1 - tj ). 


Furthermore, the field operators in the two-body term can be rearranged to 
take the form of the two-body Green's function (Eq. (4)). The double combina- 
tions of field operators already are in the form of a single-particle Green's 
function as in Eq. (3). Combining Eqs. (3) through ( 6 ) yields: 


i 



+ 



Gj (1 , 1 ’) 




T z ) G z (l,2; l’,2>) 


(7) 


= s ( 3 ) (r, 


“r { ) 6 (tj 




which is the equation of motion for a single-particle Green's function coupled 
to a two-particle Green's function. One must either approximate the two- 
particle function or write its differential equation which couples G 2 with Gi 
and G 3 . The last alternative is rarely done. Equation (7) expresses the 
manner in which single-particle propagation is coupled to and affected by the 
presence of all other electrons. 

Oftentimes it is more convenient to rewrite Eq. (7) as a purely single- 
particle Green's function equation by introducing a new integral operator, 
called the mass operator. In this new form, Eq. (7) becomes 
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( 8 ) 


a 

(l at 1 


v 

T~ 


c 

-) Gj ( 1 , 1 1 ) - jVr’ 1 Z(l,l") Gj (1 ", 1 ') = 6(1 - 1') 


The price to be paid for this simplification is that now Eq. (8) is most likely a 
non-linear equation and Z is a complicated, non-local integral operator. 

Equation (8) is particularly easy to solve in the case of homogeneous 
systems. In this case, absolute position has no significance, only relative 
separations. Thus, the Green's function and mass operator become functions 
of one variable only in the manner of G ^ ( 1 , 1 1 ) — * G . (1 - 1 1 ) . If the F ourier 
transform of Eq. (8) is taken with respect to (1 - P) for the free-field case in 
which Z (1, 1") = 0, the f ree-particle propagator is: 

G q (p, co ) = \ (9) 

CO - ± i T1 

where r| is an infinitesimal constant, specifying the contour of integration. A 
formal integral equation, known as the Dyson equation, can be built up from 
the Fourier transform of the full Eq. (8). The momentum space integral 
equation for the full interacting propagator is giverTby 


G(p, co ) = G o (p, co ) + G o (p, co ) Z(p, co ) G(p; co ), (10) 

which is formally solved as 

G (p, co ) 

G(p, <o ) = — (11) 

1 - G (p; co ) 2 (p, co ) 

The energy eigenvalues of the interacting system are ascertained by the poles 
of the full Green's function of Eq. (11). Direct substitution of Eqs. (9) in (11) 
gives the eigenvalue condition: 

^ 2 

w (p) = + Z (p, <*> ) 


which is a complicated transcendental equation caused by the dependence of Z 
on co. 


The usefulness of the Green's function method is seen in the following. 
The expectation value of any single -particle operator with respect to the 
exact interacting system is given by 



2^r / d “/ 


d 3 r 


O(r) G^r.O ; 


0 + ) 


( 12 ) 
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where the + sign signifies that the time variable in the creation operator is 
taken at an infinitesimally later time than in the annihilation operator. The 
contour in the complex energy plane is chosen such that it encloses the poles 
of all occupied states. In particular, the charge density at a point is given by 


/ 


dco (r, 0; 


In general, the type of self-energy processes represented by the mass 
operator Z being considered are those resulting from a particle of momentum 
"p propagating in the many-body field which induces a virtual density fluctua- 
tion of momentum "cf, thus reducing the particle's momentum to'p' - "q. The 
particle and density fluctuation propagates until either a new density fluctua- 
tion is excited or the original density fluctuation and particle recombine. 

These processes are best visualized in terms of Feynman diagrams as shown 
in Figure 1. In these pictures, a single straight line represents the free- 
particle propagator G (p), a double straight line represents the exact propa- 
gator G(p), the wiggly line represents a propagator for the effective inter- 
action called 


D (q) 



and the circle is called the vertex function representing non-linear complica- 
tions which result as the particle and interaction lines recombine. 
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Figure 1. -- Typical diagrams contributing to the self-energy of a 
Fermion propagating in the many-body field. 



It can be shown that the sum of all proper self-energy processes of the 
form shown in Figure 1 can be exactly summed to yield the single diagram- 
matic representation also shown which is expressible as an integral: 

r h 4 

Sfp, p°) = i / -5-^2 G(p - q) r (P.q) D(q) . (14) 

J (2tt) 4 

I p" -*q|> p f 

The restriction on the allowed values of q guarantees that the Pauli principle 
is not violated in the intermediate electron state. 

Equation (14) can be rewritten in the form: 

4 

S (plp°) = V exch (p) + i J ^4 G(p - q) p (p,q) D F (q)(v(^)) 2 (15) 

with the density fluctuation propagator identified as 
D F (q) ~ v(q) ( e (q) 1 ) . 

Equation (15) follows from the fact that a Green's function, or in the field 
theoretic language, a propagator for Poisson's equation can be identified 
which is the differential equation whose solution is the screened 1/r potential. 
In Eq. (15), the interaction propagator has been split into two parts. The first 
is a generalized exchange term; the second involves the creation and propaga- 
tion of a density fluctuation through the system. As will be seen later, the 
density fluctuations arise from exciting states within the electron hole con- 
tinuum or exciting collective plasmons. The evaluation of Eq. (14) for vari- 
ous systems will be one of the major undertakings in this thesis. 

It should also be pointed out that a Coulomb potential can be written as a 
Fourier integral: 

e I d q lq. r 

— = I —■ i v (4) e H 

J (2<r 

with the Fourier coefficient 


v(q) 


4-rre 

2 

q 


In the theory of the free electron gas, abundant usage of a dimensionless 
parameter r g prevails. This parameter is uniquely determined by the elec- 
tron gas density through the following relation: 
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No. electrons/ vol = 1 /-i- irr^ and r = r/a 


where a n is the Bohr radius. For metallic densities 1.8< r s < 5. The 


RPA is strictly valid only when r s < 1,8. Various formulas for free elec- 

tron gas quantities of interest appear below as functions of fundamental con- 
stants and as functions of r s . 


Number density: 


n 


1. 61 


3 3 

4na r 
o s 


A ' 3 


Plasma frequency: 


t0J = * (±EI 

p V m 


2\l/2 


' 12 \ 1 / Z 
“T \ Ry 


Fermi wavenumber: 


4 , 
9 " 


ITT 


a r 
o s 


3. 64 


P- = 1.72 r X 
k^ s 


A 


-1 


< E >= ! 


2 . 21 


2m 


Ry 


Thomas-Fermi screening length : 


6trne 


2\l/2 


"T- F 


2 . 22 

"772 


o 

A' 


Since many-body theory is usually tried out on materials with simple 
band structures, typically the alkalis, beryllium, and aluminum, numerical 
values for the fundamental electron gas properties are displayed below: 


20 



^~'~~'~----~-^^Material 
Pr ope 

Be 

A1 

Li 

Na 

k 

r 

s 

1. 88 

2.07 

3. 22 

3. 96 

4.87 

°_ 3 
n (A 

0.243 

0. 182 

0. 048 

0. 026 

0. 0139 

hw (R ) 

p y 7 

1.35 

1. 17 

0. 6 

0.44 

0. 321 

° -1 
k F (A l ) 

1. 94 

1. 76 

1. 13 

0. 92 

0. 747 

< E > <v 

0. 626 

0. 515 

0. 213 

0. 141 

0. 093 


1. 62 

1. 54 

1. 24 

1.11 

1.01 


At this point all the formal apparatus required to understand the following 
sections has been cited. 



Ill, FIELD THEORY OF AN INTERACTING ELECTRON GAS 
IN A PERIODIC POTENTIAL 


Generalities 

Since the investigations undertaken in this report provide basic under- 
standing of real physical phenomena occurring in real metals, it is at once 
necessary to see what effects the presence of a periodic potential have on 
many-body effects, since this is a very critical inhomogeneity. Frequent ob- 
jections arise from the fact that most many-body calculations are done using 
a plane wave basis to describe the originally non-interacting electrons. One 
then seeks the adiabatic transform of plane wave electrons when the e-e inter- 
actions are turned on. But so-called non-interacting single electron 
"eigenstates" in metals are Bloch functions and not plane waves so one should 
seek the adiabatically transformed Bloch state if one is dealing with many- 
body effects in real solids. As mentioned in section I, this problem has been 
considered to some extent by others (refs. 11-14). 

The theory to be presented in this section contains some ideas which lie 
implicit in references 11 through 14 and are here made explicit in physically 
understandable terms. The basic underlying concept used is the following. 

The system of concern has two different sorts of interaction terms in the full 
many-body Hamiltonian- - electron- electron interactions and electron-lattice 
interactions. The electron correlation part is treated first, obtaining Green's 
functions for single quasi-particles which include a self- energy term resulting 
from propagation through the polarizable field of the Fermi sea. Then the 
dressed electron is allowed to propagate through the periodic pseudo-potential 
and the self-energy of a quasi-particle in a periodic potential is obtained. 

The correlation effects are evaluated in the plane wave approximation. 
This is a reasonable assumption though for the following reasons. First, the 
procedure of Bassani, Robinson, Goodman, and Schrieffer (ref. 13) could be 
followed by expanding electron states in an orthogonalized plane wave (OPW) 
basis and then considering many-body effects on the OPW electrons. This, of 
course, makes the mathematical aspects of the many-body part of the problem 
much simpler than if an attempt to compute correlation effects on aug- 
mented plane wave (APW) electrons had been made. From a physical point of 
view, the plane wave approximation can be seen to be sensible. In the core 
region, far within the APW sphere, that region in which electrons do not have 
plane wave behavior, the electrons have high kinetic energies and thus go 
through the core very quickly. Consequently most of the electron's time is 
spent in the region in which it behaves like a plane wave and, for the most 
part, the electrons will feel correlation effects as if they were plane wave 
electrons. In the core region, the high energy electron will not have its be- 
havior disturbed much by the relatively weak correlation interaction so there 
is really little need to carry through atomic -like core effects, which are well 
treated in the Hartree-Fock approximation, in calculating many-body correla- 
tion effects; that is, many-body effects are important only where core effects 
are small so an effective separation of the two can be made. This sort of 
argument could be expressed quantitatively along the lines of the Cohen-Heine 
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cancellation theorem (ref. 48). It is worthwhile to realize that high energy- 
electrons don't correlate very much relative to the plane wave portion of the 
electron, thus permitting treatment of many-body effects before the lattice 
effects . 


The structure of this section is as follows. In the next subsection. Theo- 
retical Considerations, general field theoretic considerations are given to the 
problem of a particle moving in two fields -- the electron gas and the periodic 
potential. Expressions are obtained for the single quasi-particle self-energies 
and thus for OPW energy bands with many-body effects included. In the last 
subsection. Discussions and Conclusions, a discussion of the physical origin 
and implications of these results is given. 


Theoretical Considerations 

Consider the full electron Hamiltonian for a perfect real metal with inter- 
acting electrons: 


where 



V i i 
-T- + N V<r 


i. j = l 


'i + R j>} + X v(r i- r j> 

j i< i 


V(r^ + IC) is the periodic lattice potential and 


v(r. - r.) = 

1 J hr 


r. - r . 
1 J 



iq • ( r. - r .) 
, i J 


is the electron-electron interaction. Zero temperature is assumed and electron- 
phonon interactions are neglected. 


Consider now the following part of the full Hamiltonian: 

2 

-^- + * ^ 

i = 1 


N 2 

H„- X -3-* X 


. . r. - r . 

K j 1 i J 


N 

X 

i = 1 



iq • (r. - r.) 
1 J 


which is needed to describe the dominant many-body effects. The starting 
point is the kinetic energy operator whose eigenfunctions are plane waves. 

This allows defining of a spherical Fermi surface which is equivalent to the 
vacuum state in field theory. Since all the momentum eigenstates are filled 
beneath the Fermi surface, there can be no particle propagation in momentum 
space unless a passive electron is excited above the Fermi sphere. Then the 
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elementary excitation with momentum greater than the Fermi momentum, 
called an electron, can propagate freely, scatter off of a source, or interact 
with a field. Likewise, the positively charged elementary excitation called 
a hole, with momentum less than the Fermi momentum, can undergo the 
same sort of processes. 

Next the modification upon propagating electrons with p > p^, passive 
electrons with p < p^, and propagating holes with p < p^ which result from 
turning on electron-electron interactions should be studied. Two passive 
electrons within the Fermi sphere interact through the q'th Fourier component 
of the e-e interaction, giving two virtual electron-hole pairs. This is shown 
in Figure 2. These propagating electrons and holes can then interact with the 
polarizable Fermi sea or with the lattice potential. 

|P + ql>P F 



Figure 2. -- Interaction of two electrons beneath the Fermi 
sphere resulting in the formation of two 
electron-hole pairs. 

Consider first self-energy processes within the Fermi sea. The propa- 
gating Fermions can have their lines modified by processes such as those 
which were shown in Figure 1. The first object of interest in our calculation 
is the self-energy resulting from e-e interactions as given in Eqs. (1) and (4): 

z e ( p . p°) = i f - G e (p"- q; p° - q°) r (p, q) D(q) , 

J (2 ’ r) (16) 

It-qj > P F 

with the meaning of relevant quantities well discussed in Section II. In order 
to obtain calculable end results, approximations on the quantities appearing 
within the integral of Eq. (16) must always be made. The first of these is to 
set the vertex function T (p,q) = 1. This is one of the linearizing approxi- 
mations. It is hard to check the validity of this assumption quantitatively be- 
cause few people have been able to find manageable methods of going beyond 
this zero- order interaction. It might be hoped that since E e (p) is expected to 
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be a small correction to the single- electron spectrum, corrections to F = 1 
would be of a higher order smallness and thus not too significant. 

As an example of the type of corrections resulting from higher order 
vertex functions, consider the following. It is known from BCS theory of 
super- conductivity (ref. 49) that small corrections to the zero- order electron- 
phonon vertex function give rise to the superconducting instability; that is to 
say, vastly different qualitative changes can occur as a result of small modi- 
fications to the vertex function. Since metals which have quite normal pro- 
perties are being considered here, there are few qualitative disparities between 
intuitive theory and experimental observations. Thus it appears that the zero- 
order vertex function is a justifiable approximation since small corrections can 
lead to qualitative changes which are not required. 

With regards to the exact interaction propagator D(q), there exist various 
standard approximation schemes for obtaining dielectric functions, each with 
associated good and bad points. The Thomas- Fermi function has a simple 
functional form which is good from the computational point of view, but has no 
dynamical properties and hence describes static screening only. This is un- 
realistic. The Hartree-Fock dielectric function has the virtue of describing 
dynamic screening within the context of a very familiar model, but has the 
unrealistic property of a vanishing density of states at the Fermi surface which 
is the region in momentum space where most many-body effects occur. The 
random phase approximation (RPA) provides an exact dielectric function, but 
only for electron densities much higher than those found in real metal con- 
duction bands. The Hubbard-modified RPA function is nearly equivalent to an 
interpolation between exact high and low density dielectric functions (ref. 50) 
and thus should provide a reasonable approximation to the exact interaction 
propagator. 

To obtain the propagator for the physical electron in terms of non- 
interacting electrons, the Dyson equation is formally solved: 


G 0 ( p " 9) 

G e ( p " q) 1 - G o (p - q) S e (p - q) 


(17) 


In Eq. (17), the free particle non-interacting Green's function 


G (p - q; p - q ) = 

o r o o 

p - q - 

is used to describe electrons (holes) if 
Fermi momentum. 


1 

(I “p -"qf ) 2 / 21 * 1 ± iT i 

| "p - "q | is larger (smaller) than the 


It is useful at this point to illustrate the types of processes which are 
usually considered (at least within the RPA or Hubbard -modified RPA (ref. 7) 
to give the most divergent contributions to the self- energy loops in Figure 1. 
Some of these are shown in Figure 3. In Figure 3a, the propagating particle 
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Figure 3 — Lowest order Hartree and RPA self-energy diagrams. 



Figure 4 -- a Hubbard modification of a bubble diagram. 

T) Multiple scattering modification of a bubble. 



Figure 5 -- Graphical illustration of the dressed electron propagator 
in a periodic potential. 





interacts with a passive particle, but not strongly enough to excite it out of the 
F ermi sea, so the passive particle must go back where it started from within 
the Fermi sphere to prevent violation of the exclusion principle. Thus the 
propagating particle simply scatters off the Hartree potential of the Fermi sea 
electrons. In Figures 3b and 3c, the propagating particle interacts strongly 
enough with a passive electron to excite it above the Fermi surface thus leaving 
behind a lower momentum hole. The pairs propagate and eventually annihilate 
by creating a new pair or by interacting with the initial pair. 

Here the relevance of the Hubbard correction to the RPA with regard to the 
periodic lattice problem should be noted. Hubbard has suggested that the RPA, 
an inherently high density approximation, comes closer to the truth for metallic 
densities if one allows the propagating electron and hole forming the pair to 
interact; that is, the motions of the electron and hole are correlated. He thus 
replaces uncorrelated bubble diagrams by the diagram shown in Figure 4a in 
which the electron and hole interact once through a Thomas-Fermi screened 
interaction. Recently Crowell, Anderson, and Ritchie (ref. 51) have investi- 
gated this problem in relation to positron annihilation theory and have found 
that relatively good agreement with experiment could be obtained if instead of 
a single screened interaction, a multiple scattering ladder series was used to 
describe the electron-hole attraction. It is significant that the series con- 
verged extremely rapidly, requiring only three or four multiple scatterings 
within a single bubble. 

In the case of the periodic lattice, another effect is present which may 
tend to cancel the electron-hole attraction. The lattice core potential is 
attractive for the electron and repulsive for the hole. Thus once a pair is 
formed via the e-e interactions, the electron and hole will feel an effective 
repulsion as a result of the core. At worst this will minimize the need for a 
Hubbard-type correction. Although no calculations as to the quantitative im- 
portance of this effect are offered, some conjectures may be made. It is 
indeed possible that, under certain conditions, a balance between the effective 
lattice induced e-h repulsion and the direct e-h attraction could be obtained. 
Under these circumstances the electron-hole pair would have quite a long 
lifetime and could be recognized as some sort of many-body excitation re- 
sembling a virtual exciton. One manifestation of such an entity might be an 
effective dipole potential around the ion cores giving rise to a mild discon- 
tinuity in the usual muffin tins. For the time being this is left as an interesting 
problem to contemplate and as a partial argument for using uncorrelated 
bubbles in real metals. 

Next one might ask how the physical Fermion lines appearing in a diagram 
describing a many-body process would be modified if these quasi-particles were 
put into the field of a weak periodic pseudo-potential. 
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xhe periodic pseudo-potential can be written as 
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( 18 ) 


where W(|q|) is the spherically symmetric OPW form factor, 

S(q) = 26 (q" + K) is the lattice structure factor, and K is a reciprocal 
K 

vector (ref. 52). In other words, the Fourier transform of the electron- 
lattice interaction is given as 

V® = W ( fq"| ) S (qj. (19) 


A quasi-particle propagating in this field will not have a true self-energy 
because the lattice potential plays the role of a distribution of sources and 
sinks rather than a polarizable medium through which the particle propagates. 
However, one can still define a mathematical object which serves the same 
function as a self-energy operator. This object expresses the way in which 
the propagation of a quasi-particle is modified by the static lattice "field". 

In Figure 5 are shown diagrams corresponding to the simplest interactions. 

If the particle lines interacting with the lattice field were those for bare 
electrons, then consideration of the two events shown in Figure 5 would be 
sufficient to generate OPW bands (ref. 52). The first figure describes the 
effects of the q = 0 component of the OPW form factor which has been canceled 
by the q = 0 component of the electron- electron interaction. The second event 
is a quasi-particle of momentum p - q scattering into an intermediate, un- 
occupied state ’p’-’q+^T , and then returning to its original state in a manner 
quite similar to the self- energy process shown in Figure 1. An effective lattice 
self- energy or level shift operator is then identified by 


2^ (p ~ q; 



V(-’q') G e (p -”q + q';p° - q°) V (q 7 ) 


I T - q'+'q' I > P F ( 2 °) 

in analogy with Eq. (16) where G e (p - q + "q ' ; p° - q°) is the propagating 
physical electron Green's function given in Eq. (17). Substituting the periodic 
potential given by Eq. (19) into Eq. (20) and performing the integrals yields 
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2L (P - D 


Z I W(K)| 2 G (p- K; p° - q°) 

K 

1 "p - <T + K 1 > P F 


( 21 ) 


Furthermore, considerations in this report are limited to points in the first 
Brillouin zone away from edges and corners so that, at most, only two OPW's 
need be considered. 


The formal procedure of writing and solving a multiple scattering t ma- 
trix or ladder sum for a quasi-particle Green's function modified by the ef- 
fective lattice self-energy Z T can be used to obtain 

JLi 


(p - q) 


Ge (p - q) 

1 - G e (p - q) Z L (p - q) 


( 22 ) 


Using Eqs. (17) and (21) in the above yields: 

g l (p - q) = ^[(p 0 - q° - E 0 (p"-"q) - £ e (p’-'q; p° - q° )) 

- | W(K)| V(p° - q° - E o (^-T+K) (23) 

- S e (^--q + K; p° - q° ))] - 

If r (p, q) is set = 1 as discussed earlier, this Green's function describing a 
quasi-particle propagating in a periodic potential is placed into the original 
self- energy expression, a slight amount of rearrangement is done, the elec- 
tron self-energy, valid for all values of p, is: 


: Tot^ P ’ P * " 


it * d 4 q o. 

1 / — "4 D (q; q ) 
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IT -"q| > p F 
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£(p° - q° - e q (p -"q) - s (p - - q; p° - q°) ) 
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x (p° - q° - E o (p - q + K) - S (p - q + K; p° - q°) ) 


- |W(K) | 2 ] 


(24) 
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The positions of the poles of the integral are determined by the solution of the 
secular determinant: 

p° - q° - E q (p'-'q) - Re 2 e (p - q; p° - q°) W(K) 

W(K) p° - q° - E Q (p- - -q + K) - Re S (p - - q + K; p° - q°) 

with ReZ e the real part of the polarization self- energy operator. It is known 
that E e is very nearly constant with variations in momentum .(ref . 9). Thus: 

2 e (p - q; p - q ) = 2 e (p - q + K; p - q ) 

It can be shown that this assumption implies that there is a correlation hole 
around each electron, within which the charge density variation has the usual 
jj(k c r) spherical Bessel function form. Equating the self-energies yields: 

-q° = J (E o (p - q) + E o (p - -q + K) - 2p° + 2 Re Z (p - q) ) 

±(5 (E o (p"- J q) - E o (p--q + K) + ^(p - q; p° - q°) (25) 

-S e (P - q + K; p° - q\) Z + |W(K)| 2 ) 1/2 

- - a 

To be completely rigorous, the cut representing the continuous distri- 
bution of poles of the interaction propagator D(q) should also be considered. 
However, in a spectral decomposition, the two poles for q° will be quite sharp 
and not obscured by the background. The effects of the D(q) cut can be ab- 
sorbed in an effective renormalization constant, called Z ? , which usually sat- 
isfies 0.75;$ Z 2 $ 1 in electron gas problems (ref. 9). This renormalization 
decreases the magnitude of the residues at the two principal poles since some 
of the spectral weight now appears in the cut. 

To perform the q° integral of Eq. (24), the contour shown in Figure 6 
was chosen. The pole corresponding to the state in the second zone is marked 
with a broken x. If work is done in an extended zone scheme, this pole is on 
another Riemann sheet. Since in the reduced zone scheme all poles are on 
the same sheet, the higher zone poles must be included by introducing band 
indices. The q° part of the integral is of the form: 
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(26) 


where an = q° is the solution for intermediate states in the n'th band of the 
hopelessly non-linear Eq. (25). Referring to the denominator of Eq. (24) 
and setting the E's equal gives: 


g(q°) = 


2^ (p - q) - Z. (p - q) (2 p° - 2 q° - E (p“ - "q) 


o 

o o x 


- E o (p - q + K) ) + (P° - q°) 2 - <P U - q”) (E Q (p - q + K) 
+ E (p"-"q) +E n (p-q) E (p - q + K) - | W(K) | 


(27) 


C onsequently: 
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and one of the Ward (ref. 53) identities is used: 


A (- q°) = 1 


^ z (- q ) 

8 q° 


then Eq. (29) takes the form: 

(p - q; p° - q°) 


= l - A (q°) - p° (i - A (-q°) ) . 

8 q 


8 q 


( 30 ) 


Taking the vertex function /\ (q ) = 1> the usual first approximate which is 
equivalent to linearization: 


g'(q°) - (-( E 0 (p - q) + e g (p -q" + k) - 2 ( P ° - q°) 

+ 2s e (p - q; p° - q°)) . (31) 

Now, in a similar manner, if the expansion, Ward's identity, and the unit 
vertex function for treating difficult E 's in Eq. (25) are used, then 
gi (qO ~ a n ) becomes simply: 

g'(a n ) = ± 2 ^( E 0 (P " q) - E 0 (p - ^q + K)) 2 + |W(K)| 2 y /2 (32) 

Thus through the complex integration of Eq. (26) the total self-energy 
operator, Eq. (24), becomes 


Srot ( P ; P 0) = -Z Z 2 S / a n )(p° - a n - E Q (p - q + K) 

n=l, 2 J (2tt) 


| p - q‘| >P F 
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(33) 
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Putting Eq. (25) for into the numerator gives the result: 

^Tot ( p;p ) = " 2 Z zJ ^^3 [ D ^ q; a l^ " D ^ q; °?) J 
|p‘-"q| > p f- 


-i Z zf 7^3 [ D <* a i> + D <* ° 2 > 
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|P~ - "q | >^ P F 
Jp - - ■q + K| > p F 


(34) 



E Q (p - "q + K) - E q (p - q) 
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J <E 0 (p - q) - E o (p - f + K>) 2 + 

W(K) 
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1/2 


Note that the only place in which p appears in the right hand side of Eq. (34) 
is through Eq. (25) determining the value of a n used in the interaction prop- 
agator D(q~; a ). The energy spectrum is given by 


E (P) = e opw <P> + ^Tot <P : E) 

with the first iteration 


E (P) = E opw (p) + S rot (p; §n> ). (35) 

These are amazingly simple and quite understandable results considering the 
vast complexity of the system investigated. The meaning, interpretation, and 
discussion of the present results appear in the following section. 


Discussion and Conclusions 

Through a rather involved mathematical excursion, a very simple 
general expression for the self-energy of an electron in the many-body po- 
tential of an interacting electron gas, plus a periodic lattice potential, has 
evolved. The result is quite general, subject to the folliwing condi- 
tions. Since the RPA has yet to be invoked, this not necessarily a high 
density theory. To perform numerical calculations would require an explicit 
expression for the interaction propagator D(q) = v("q) /e (q). If plane wave 
many-body theory can give a good e (q), then Eqs. (34) and (35) can give a 
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correspondingly good self-energy if there exists a periodic potential, at least 
for materials and points in the Brillouin zone which have been previously 
treated by two non-interacting OPW electrons. First attempts will be made 
by using the RPA or Hubbard- corrected RPA to see how significant many- 
body interactions really are (ref. 54). It would also be interesting to split 
up the density fluctuation propagator into a plasmon and pair propagator 
as was done by Hedin et al (ref. 10) and see how the lattice affects the pla_smon 
part of the self-energy. 

A look at Eq. (34) prompts the following comments as to the structure of 
this end result. The first simple integral is just one-half the self- energy for 
an electron gas with no periodic lattice except for two modifications. First, 
the interaction propagator is taken at the energy a n of Eq. (25) rather than 
at the free electron energy because the lattice shifts the free electron en- 
ergies. Second, the fact that work is being done in the reduced zone scheme 
is made manifest by admixture with the OL 2 pole. The second integral pro- 
vides the dominant lattice contribution. In the empty lattice test where 
W(K) =0,0;= the free electron energy and the second integral becomes one- 
half the free electron self-energy with the modification from the admixture 
with the Qf2 P°l e * However, this contribution exactly cancels the term in 
the first integral so the required limiting expression is precisely obtained. 

That Eq. (34) gives the required value in the limit adds confidence to its 
overall validity. 

The origin of the lattice contribution can be explained as follows. Initially 
an unperturbed Fermi sphere exists. Two passive electrons interact as 
shown in Figure 2 and are excited above the Fermi sphere. If they have 
interacted through a large enough momentum transfer q, then one of the vir- 
tual electrons may have enough momentum to satisfy the Bragg condition 
2(|T + "q) • K = K . If one of the virtual physical electrons (i. e. , corrected 
for many-body effects) is Bragg-reflected, then the remaining hole will not 
have any electron to annihilate with and thus will have a longer lifetime to 
make its influence felt. The same is true for the reflected electron. 

Consequently both the electron and hole will remain in high energy states 
for longer times than are advantageous for reduction of the total energy. Thus 
the dominant important lattice effect of Bragg reflections serves to decrease 
the absolute value of the self-energy. Since the self-energy is a negative quan- 
tity, the lattice increases the exact total energy of the system relative to the 
free electron total energy when evaluated at a^. 

It is worthwhile commenting upon the procedure to be used for obtaining 
numerical results based on the theory presented herein. To calculate many- 
body OPW bands, pick a particular symmetry line and point on this line. 

This determines the value of p used to calculate a n given in Eq. (25). Know- 
ing the symmetry line of the particular lattice considered, one can then ob- 
tain tabulated values for the OPW form factor W(K)(ref. 52). Using the numer- 
ical tabulations of Hedin, one can also obtain values for required in 
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Eq. (25). Expressions for DCqja') exist in the literature within the various 
approximation schemes for the dielectric function as previously discussed 
(ref. 6). Thus everything needed to perform numerical integration of 
Eq. (34) exists. Sample calculations to see exactly the lattice effects on 
the many-body theory* are in order. 

Calculations of some experimentally relevant quantities such as 
electron specific heats, density of states at the Fermi surface, and band- 
widths, or equivalently, effective masses are also in order. It will be 
interesting to see what percentage of the electron effective mass, usually 
attributed to single-electron interactions with the periodic potential, 
actually arise from many-body interactions in a periodic potential. 

An important point to be picked up from this work is the following. 
Electronic behavior in solids is determined in large by three distinctly 
different phenomena which can be treated separately. The first of these 
is concerned with high kinetic energies of electrons in the core region. 

One can project out of the true wavefunctions the wiggly portion which is 
well treated in an atomic-like Hartree-Fock approximation. The second 
point is concerned with the possibilities of Bragg reflections at zone faces. 
This is dealt with by taking linear combinations of the wavefunction which 
remain after projecting out the high-energy components. The last point is 
related to the many-body effects which, for the most part, affect the low- 
energy or plane -wave components of the exact wavefunction. By doing 
each problem separately and as exactly as possible and then combining the 
results in a self-consistent manner, one can finally arrive at potentially 
useful results which should provide a better understanding of real electronic 
behavior in real metals. 

Thus, the first plateau in this study of inhomogeneous electron gases 
has been attained. An expression has been obtained for the momentum- 
dependent self-energy of an electron in an interacting electron gas with a 
periodic pseudo-potential. Although the importance of this correction to 
single electron energy band theory may be of significance, as will be seen 
in the near future, it appears that the lattice offers no quantitatively new 
feature which needs to be included for present purposes. It is thus con- 
cluded that lattice effects may be of significance in the many-body correc- 
tions to energy band theory, a quite precise theory, but are not of first 
order importance in the inhomogeneous electron gas effects to be treated 
in the following sections. 


*P. deCicco (personal communication) 
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IV. DIELECTRIC FORMULATION IN AN ARBITRARY, 
NON-INTERACTING ELECTRON BASIS 


Generalities 

As has been stressed previously, the central rallying point of many- 
body effects can be considered to be the time-dependent linear response func- 
tion of the electron gas to a driving force. This response function is usually 
presented in the form of a dielectric function, susceptibility, or polariza- 
bility. The simplest derivation of the RPA dielectric function for a homo- 
geneous electron gas has been through the self-consistent field theory of 
Ehrenreich and Cohen (ref. 8), further elucidated by Brout and Carruthers 
(ref. 54). 

A general derivation of the RPA dielectric function for an inhomogeneous 
electron gas, within which the originally non-interacting electron wave func- 
tions are expandable in a plane wave basis, is now presented. In the limit in 
which each wavefunction is the "sum" of only one plane wave, the homogen- 
eous case, the derived dielectric function reduces to the ordinary RPA re- 
sult. With the aid of a new image technique, the general inhomogeneous 
electron gas dielectric function takes on a simple form in the surface region 
also. With this imaging technique, the inverse dielectric function appropriate 
in the surface region displays a pole at the surface plasma frequency. 

Theoretical Considerations 


Following EC, the Liouville equation describing the time evolution of the 
single particle density matrix is considered: 


. 3 


o_2 
3 t 


[h,p]. 


(36) 


The Hamiltonian for the system is written as 

H = H + V(r,t) 
o 


(37) 


with 


H = KE + V (r) 
o o v ' 


defining the complete set of non-interacting electron eigenfunctions | ip \ 
which are assumed to be plane wave expandable. The exact eigenfunc- 
tions satisfy 
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Nk) 


= E(k) i r 


l*tr > 


(38) 
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with 



and T the coefficient matrix. The self-consistent potential V(r, t) is 
determined by the applied field plus the possibly time -dependent induced 
screening charge. The operator p, representing the single-particle den- 
sity matrix, is expanded as 


P = P 0 + Pi 


(39) 


such that 


Pol k) = f o 1 ^ 

where f is the Fermi function, 
o 

Considering that the non-interacting Liouville equation is 


(40) 
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3 p 

F o 

3 t 



(41) 


and linearizing by neglecting 
[ VCr,t), p ] 


the term 


in the full commutator of Eq. (36), allows determination of the equation of 
motion for pj , the measure of departure from the equilibrium distribution. 
Substitution of Eqs. (37), (39)> and (41) into the linearized Eq. (36) yields: 


3 p x 
1 “3T 



[ p o . 


(42) 


Taking matrix elements of Eq. (42) between the time -independent states 



and utilizing Eqs. (38) and (40), then yields: 

<^klpl l^r + q ]> = < E < r > - E ' r+l J» ^r'pl^k+t ^ 

+ ( f o < E r + q ) - f o < E k» )\ * r i v < r -‘> 1 * ir+ £>• (43) 

In keeping with the sign convention of EC, the self- con si stent potential 
is written as 

V(7,t) = £ V(q,t) e" iq * r . 

and is related to the charge density through Poisson's equation: 

„ 2 __ . 2 

V V = -47re n 

or in momentum space. 


4ire 

V(q,t) = — 2- n- 

n ~ 


(44) 


with 


n < r > = Z n ^ e_iq ' r * 


(45) 


The expectation value of a single particle operator is given by 

<V> = Tr ( p O) 

where p is the density matrix and Tr is the operation of taking the trace. 
In the present case, the density increment operator is of interest. 

Pi = 5 <* op - J?) - 

The trace is written as the sum of diagonal elements of the matrix of the 
operator between the exact single-particle states j . Following this 

procedure, the density increment is: 


n(T) = Z < ( 4 , k' i 6(? o P - T > pi i * k>> • 

vi \ ' 


(46) 
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Now inserting the exact eigenfunction representation of the unit operator 

i = I 1 <■** 1 

k 

into Eq. (46) yields: 

nM = I 6 <Kp | 4> k I Pi I ^k-> * (47) 

k.k' 

If a complete set of eigenfunctions of the position operator is introduced: 


J d 3 y | y > < 


1 = 

< 

then Eq. (47) becomes: 
n(r) = 


) = J d3 y ^ 0 £7 | y (y - y I < ^l p i 

k,k ' (48) 


with the transformation functions 


r k,k"' e 


ik"',y 

■ ™ P J 


<^k' \y> = ? 


(49) 


Performing the y integration and using the transformation functions allows 
Eq. (48) to take the form: 


T\ - Z r- 2 . - r— i(k"' - k").T , x 

__ r k',k” r k,k'" e <^klpi l^k,> * 


k, k' 
kj' "k' n 


(50) 
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If one calls k' = k + “q and replaces the sum on k' by a sum on q, then: 

n(r) = Z e- i? - T Z 

■q "k, k"' ’ 

k m 

<*klPl I*ir + q>- < 51 

In order that Eq. (45) and (51) be compatible, it is required that 

tT' " - IT' 1 + q = 0. Thus the Fourier coefficient of the density increment is 


-P* i(k~'" - k" + q) . r 

1 IT + "q,k" k,k"‘ 


n(q) - X r tT+qf,k i ' r 'k,k"-'q ( I P 1 l^+q“ ) . (52) 

k.k 1 ' 

Returning to Eq. (43), it if it assumed that the time dependence of 
goes as 

-icot 

Pj - e 


and if the excitation frequency is defined as 
E(k~ + "q) - E(k) = co-~(k), 

then 

(f (Es-,—) - f (EH) 

v o v k + q o k' ; 
oj + u>-?-(k) 



- ( I V < r> I * k + q ) . 

The matrix element of the self-consistent potential is written as 


(53) 
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Since the plane wave matrix element is 6(k - q' - k ), Eq. (54) becomes: 


(** l V( r) I ^ k + q) Z rl r_ _ _ v(q'). 


‘k^.'J' k ’ k V k + q ,"k * v + q ’ 


(55) 


Combining Eqs. (52), (53), and (55), and changing variables k" — ► k 1 and 
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(56) 


Classical electrodynamics is now used to obtain an expression relating 
the dielectric function with the polarization charge. The polarization is re- 
lated to the induced charge through 

V * P = en 

or in momentum space 

-iq * P-»- = en-*-, 

q q 

The polarization is also given by 

p_ = «(qA) ~ 1 

q 4 t r q 

with 


eE a~ = -iq 
HL 

Solving for the dielectric function: 

*(q,co) = 1 + v(q) — . 

V(q,t) 


(57) 
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If Eq. (56) is linearized, as in the spirit of the EC theory, by including only 
the q' =. q term, then Eq. (57) is: 


~ f (E- -) - f (E-r) 
« (*T, w) = 1 + v(q) X — 2 ll_ 


(58) 


x J X r ^ r 

j k + q , k 1 k 5 k ' _ q k > k " 

k\k" 


k + q, k " + q 


This is a completely general result for the linear dielectric function of any 
inhomogeneous electron gas, in which the non-interacting electron wave- 
functions are expandable in a plane wave basis. If the bracketed term were 
unity, then Eq. (58) would be the precise expression for the RPA dielectric 
function in a uniform electron gas. The utility of Eq. (58) depends upon the 
simplicity of the original non-interacting basis. Clearly if the F matrix is 
complicated, requiring the superposition of many plane waves, then the re- 
sulting e , if calculable, will be very complicated. 

In the following section, it will be shown that Eq. (58), together with an 
image technique, can produce a relatively simple expression for the RPA 
dielectric function in an electron gas with a surface. 

Surface Dielectric Function 


One of the more interesting and experimentally relevant examples of an 
inhomogeneous electron gas occurs in the surface region of a metal. As is 
always the case in these sorts of problems, the destruction of translational 
invariance reduces the symmetry of the problem to the point where any 
quantitative theory can be attempted on highly idealized models which 
accentuate only the most dominant features of the physical system. 


The problem of surface properties of electron gases has always been 
handled within the context of the following model (refs. 33-35, 39). A po- 
tential well with a perfectly reflecting infinite surface barrier perpendicular 
to the z direction is assumed. This allows one to work with the following 
non-interacting wave -functions: 


ip (r) 



^ (r) = 0 


Z 5: 0 


(59) 
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In general, people like to think of decaying exponential wave-functions for 
z ^ 0. However, if one considers the nodal plane z = 0 as a mathematical 
surface somewhere outside the physical surface, then the charge density 
given by the Eq. (59) wavefunctions 


n 
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4 

(2t r) 3 



. 2 
sm 


k z d 3 k 
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n. , 
int 


sin (2ppr cos 0) - (2p^r cos 9) cos(2ppr cos O)'' 1 

1 ‘ 3 

(2ppr cos 9) 


(60) 


is not significantly different from that obtained with decaying exponentials 
and the physical surface at z = £, £ > 0. As Bardeen has suggested, 

% - 0(1 A). 

The coordinate system used in Eq. (60) is illustrated in Figure 7. In 
Figure 8, Eq. (60) is drawn for various radial lines into the electron gas. 
Note that if 


Ppr cos 0 

was taken as the ordinate, then a single universal curve would be applicable, 
that one for 0 - 0 in Figure 8. 

To proceed an expression for the dielectric function of this model is 
needed. The coefficient matrix F is obtained by taking the Fourier trans- 
form of the Eq. (59) wavefunctions. Thus: 
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Figure 7. -- Coordinate system used to describe 
electron gas surface effects. 



tion of position near an electron gas 
surface. 
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Figure 9a. -- Model for describing the lattice re- 
sponse of a surface atom in terms 
of the response of an atom in an in- 
finite area. 




Figure 9b. -- Model for describing the dielectric 
response of an electron gas with a 
surface in terms of the response of 
an infinite electron gas. 


which indeed is more complicated than one desires for use in Eq. (58). 

Thus the motivation for finding a simpler technique definitely exists. 

Goodman (ref. 55) has devised an image technique for describing the 
lattice dynamics of surface atoms. How this approach can be modified to 
describe electronic surface properties will now be described. Briefly 
what Goodman does is the following. He realizes that a surface atom in 
the semi-infinite lattice responds to a point impulsive displacement from 
its equilibrium position in a manner which can be related to the infinite 
lattice response. In Figure 9a the situation is illustrated. The figure 
shows that the response of the semi-infinite lattice to an initial surface 
atom displacement is exactly the same as the response of the infinite 
lattice to an equivalent displacement of two atoms on either side of the 
mathematical surface such that initially there is no net force across the 
surface plane. It is clearly much easier to calculate response functions 
for the infinite lattice since one does not get into the problems of com- 
plicated Fourier transforms as in Eq. (61). 

Goodman's thinking is also applicable to the surface electron gas 
response or dielectric function. Consider Figure 9b which displays the 
electron gas analog. This picture suggests that the response of a semi- 
infinite electron gas to some perturbation a distance z ' into the gas is 
exactly that of an infinite gas with the same perturbations at z 1 and -z 1 
such that no electric field lines cross the mathematical surface. Although 
an example in which the charge density has a node at z = 0 has been drawn, 
this is not a necessary condition. The image technique is applicable 
wherever the mathematical surface is defined. The charge density will 
always be continuous at z =0, but only in a few cases will the derivatives 
be so. This makes for complicated Fourier transforms or F matrices, 
but does not vitiate the theoretical validity of the model. 


The beauty of the infinite gas with a nodal plane at z =0, as shown in 
the right hand side of Figure 9b is that the sine wavefunctions 
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valid for all space, are easily built up from a small sum of plane waves, 
two to be exact. The T coefficients needed in Eq. (58) take the simple form 
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where, as usual, the subscript t signifies vector momentum transverse to 
the surface. It is important to realize that Kronecker delta functions are 
used and always implied in future sections since, at times, sloppiness in 
the mathematical treatment of the Dirac delta function limit becomes 
evident. When these coefficients are used and the Tc' and ET n sums in 
Eq. (58) are performed, the bracketed term becomes: 
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Equation (64) lends itself to physical interpretation. All the 
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terms just eliminate the infinite wavelength component of the electron gas, 
a necessary requirement if one imposes the constraint of a nodal plane some- 
where. This restriction adds mathematical rigor to the present model but 
offers no computational problem. Similarly the third delta function requires 
one to avoid the point k - - q z when integrating on k. Since this point, for 
any given "q, is non- singular , one can formally perform the k integration as 
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(65) 
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This is valid since f(k) is smoothly and not piece-wise continuous at 

k = -q and non- singular. Thus Eq. (62) is simply 
z z 

The q z = 0 delta function introduces a coupling between the spherically 
symmetric dielectric function and the surface boundary condition as is 
better illustrated by writing the dielectric function in its final form: 
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which, by some simple manipulations involving a different way of expressing 
the Pauli principle constraints, can also be written as 


e (q, co ) = 1 
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introducing the notation for the polarizability. 

The zeroes of e (q", co ) ascertain the eigenvalues of the electron-hole pair 
continuum of excitations and also the collective plasmon state. In the long 
wavelength limit, as q - 0, it is easy to show that Eq. (67) gives the plasma 
frequency since 



| k + q |> p F 
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As q — 0, 
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in the denominator. In the numerator, the q * k term gives zero when 
integrated over k. Thus, going to an integral, the surviving term is 
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With v(q) = 4ue / q , the k integral giving the number density. 
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which is the standard volume plasma frequency. 


(69) 


Now a method for transcribing these results, by use of the image 
technique, into a form suitable for surface calculations will be presented. 
As was shown in Figure 9b, the response of the semi-infinite electron 
gas to some driving force is equivalent to the response of the infinite gas 
to the same driving force and its image. On a microscopic scale, con- 
sidering the driving force to be the e-e interactions. Figure 9c shows the 
physical situation. The response of an electron at a' to an electron at b 1 
in the surface problem is mathematically equivalent to the response of 
an electron a to electron b and its image c in the infinite sine wave elec- 
tron gas. Thus, if one wishes to know the response of electron a to elec- 
tron b at the surface where z = 0, then the imaging technique says that 
one wants the response of a to 2b, that is, a responds as if the surface 
electron charge were twice the standard electron charge. Conversely, 
if one wants to know how the surface electrons respond to all other elec- 
trons in the bulk, then this situation is the equivalent of a surface elec- 
tron of charge e responding to interior electrons of charge e/2. If this 
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Figure 9c. - -Example of equivalent force fields in the semi-infinite 

and infinite electron gas. 

modification is introduced into Eq. (69) for the plasma frequency, that is. 



then the pole for the collective oscillations of the surface electrons appears 
at 


W P1 



a result which has been known for a long time, but not derived in this manner* 
(refs. 36-39, 56). The main reason for demonstrating the surface plasmon 

now is to show that the theory developed so far is consistent with already 
existing theory. Consideration of the surface plasmon in relation to correla- 
tion effects upon the surface barrier will be covered in section VIII. 

It should also be noted that e in the surface problem is really quite Z- 
dependent. This fact has been hidden in the sense that we are looking only 
at the z = 0 case which simplified nicely. The problem of z dependences 
is contained in the position of electrons b and c in Figure 9c. However, even 
in the general case, the z dependence will not be drastic. The poles cor- 
responding to the pair continuum will be shifted an amount less than the level 
spacing, a negligible effect. The principal influence will be felt in the collec- 
tive oscillations. As one goes into the metal, the bulk plasma oscillations 
will begin to appear. There will formally still exist two poles at the surface 
and bulk plasma frequency, but the residue at the surface plasmon pole rap- 
idly decreases as the metal is entered and consequently the residue of the 
volume pole rapidly rises from zero at z = 0 to unity. 
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One final specific example, which will be the topic of section VI, should 
be mentioned. The problem of static dielectric screening of a point impurity- 
in the surface region is ideally suited for the theory just developed. A glance 
at Figure 9b shows that the screening of an impurity of charge e at z = 0 in 
the semi-infinite electron gas is exactly the same as the screening of an 
impurity of charge 2e at z = 0 in the full electron gas, an amazingly simple 
but correct outlook on a potentially complicated problem when viewed from 
too hasty a position. 


Now specific calculations using this formalism may be carried on. In 
the next section, the notions of impurity screening are developed, this being 
the simplest case of an inhomogeneous electron gas. Then in section VI, the 
surface dielectric formalism of this section and the impurity formalism of 
section V are coupled such that a surface impurity theory is produced. Finally 
in section VII, a many-body theory is developed in terms of the surface dielec- 
tric function, which yields surface many-body corrections to the surface po- 
tential theory of Loucks and Cutler (ref. 3 5). 
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V. IMPURITY SCREENING THEORY 


Early in the lifetime of modern many -body theory, the problem of static 
RPA impurity screening was treated by Langer and Vosko, LV (ref. 19). In 
this section, a Green's function formalism, based upon the Thomas -Fermi 
atom theories of Baraff and Borowitz, BB (ref. 16) and DuBois and Kivelson, 

DK (ref. 17), is developed. It is seen that to obtain the standard accepted 
results, certain approximations must be made which are never discussed in 
the literature. The formalism established in this section is then used in the 
following sections to describe surface impurity screening effects. 

Basically, the area of interest concerns a point positively charged impurity 
atom or ion immersed in an electron gas. This impurity then polarizes the 
gas as the free electrons move about to screen the effects of the impurity. A 
cloud of polarization charge is thus built up around the impurity. The present 
aim is to obtain an expression for the induced charge density as a function of 
distance from the impurity site. This is the simplest inhomogeneous electron 
gas problem, the inhomogeneity arising from the fact that the static impurity 
destroys the translational invariance. Thus the role of absolute as well as 
relative electron position becomes important. 

The starting point is the equation of motion for the single particle Green's 
function, given by Eq. (8) and rewritten here: 

^8T' +h(1) ) G 1 (1 - 1 ' ) - f dZ £ (1. 2 > G 1 < 2 > 1 ■> = Ml - 1') (70 ) 

where h (1) is the kinetic energy operator. 

Equation (70) may be expressed as a pair of equations 

/k( 1, 2) d2 G x (2, 1') = 6 (1 - 1') (71) 


and 

K(1 ’ 2 > +h<1) 

In these equations, the mass operator E contains all the information about e-e 
and electron-impurity interactions. The next step is to Fourier transform. It 
is here that one sees the difficulties which arise as a result of the impurity 
field. If one transforms Eq. (70) with respect to relative coordinates 


6 (1 - 2) + i £ (1, 2). 


(72) 
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r l2 ~ r l * r 2 


and the time interval t^ - keeping the remaining 


R = 


r l + r 2 


fixed, then Eq. (71) becomes 

J K(l, 2) d2 G 1 (2, 1') e* 1 p,r 12 e lai(tl ' t2) d 3 r 12 d(t 1 - t 2 ) = 1. 

( 73) 

It is desirable to express this integral in terms of the individual Fourier 
transforms of K and Gj with respect to 1 - 1'. This is a problem since K 
itself is position-dependent. However BB and DBK have provided a formal 
method for doing this by working in a mixed representation, that is, a 
representation in which K and G are partially momentum - and partially 
spatially- dependent. The mixed representation transforms of interest 
are defined as 


K (R,P) - y*d 3 r l2 dt l2 K (1, 2) e 


-ip. r 12 e iwt 12 


and 


G (R,P) =J d 3 r l2 dt j2 G (1,2) e~ ip * r 12 e lutl2 


(74) 


(75) 


The K operator can be evaluated by substitution of Eq. (72) in Eq. (74), pro- 
ducing: 


K (R, 


L,P) = JV r 12 dt 12 [iS -hfR+irj] 6 (r^) 6(t 12 ) e 
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To go further, one must consider the structure of the mass operator. 
Sham and Kohn (ref. 14) have introduced a mass operator for a possibly 
non-uniform electron gas which splits up the inhomogeneity effects into 
an average local potential and a short-range, non-local potential reflecting 
the detailed structure of the inhomogeneous gas: 

2(R + }r 12 , E) = <MR + ir u ) 6 (? 12 ) 

+ M(R + ~ T 12 , R (78) 

- i-r 12 ; E - 4 (R) ). 

In this expression, 4> is the electrostatic potential at the average position 
R which an electron feels as it moves from 'r. to "r^. cj)(R) does not depend 
on the details of how <j) (R + r) varies within tne short interval r~^ to "r^. M 
is a short range function which tells how an electron with the shifted energy 
E - * (R), the average energy in the interval r\ to"*r 2 , picks up the details 
of the density variation. Stated another way, it tells the chances that, if an 
electron is known to initially be at 

R ~ 2 r l2’ 


it will end up at 



12 ” 


It is a measure of the correlated motion of an electron through an envrionment 
of non-uniform density electrons. If this function were an instantaneous spatial 
average of two particle correlations between particles at 


R + 


2 


r 


12 


and 


then in the Hartree approximation, M would be a constant. Work will be done 
within this pseudo-Hartree factorization of single-particle, self -cor relation 
effects. Essentially this approximation is equivalent to the linearizing as- 
sumptions which LV used to do a screening calculation. If the induced charge 
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density is much smaller than the initial unperturbed density, then the in- 
homogeneity effects, as measured by M in Eq. (78), should be very small 
and the linear or pseudo-Hartree approximation thus valid. 

Inserting Eq. (78) in (77) and doing the r, ? integral with M = constant 
yields: 

K(*R,"p) = E - h(p) - 4>(R) (79) 


which sort of looks like an inverse mixed representation Green's function. 


To relate K with the physically relevant Green's function G, BB have 
shown that there is a differential operator 9, defined below, that connects 
G and K (ref. 16). Their result is stated as follows: 

9 [k (R, 5*)), G (R, ^0))] 


= lim 


R 




P' 


“*■ P 


j.K(R,p W )G(R*,^ u ) - 1. 

(80) 


As has been shown by BB and DBK, it is difficult to use the 9 operator in 
anything beyond the zero- order term in the formal expansion of the ex- 
ponential operator. It should also be pointed out that the method of gradient 
expansion of the charge density for doing inhomogeneous electron gas prob- 
lems, introduced in a somewhat phenomenological sense by Kohn et al (refs. 
29-31), would have the preceding line of reasoning as its formal justifi- 
cation. Also worthy of note is the fact that treating inhomogeneity correc- 
tions in the manner of gradient expansions will not produce Friedel os- 
cillations in the induced electron density which screens the field producing 
the density gradient. This is why the BB and BKD corrections to the Thomas- 
Fermi atom cannot, even in a formal sense, reproduce shell structure of 
the atom. Kohn and Sham (ref. 39) have made an effort to get around this 
problem by using WKB wavefunctions which necessitate a classical turning 
point and thus partial reflection at points in the inhomogeneous medium. 
Consequently interferences are set up in the inhomogeneous gas which is 
screening the field producing the inhomogeneity. The direct result of this 
phenomenon is the presence of Friedel oscillations. 

Returning to Eq. (80), one notes that to zero order in homogeneity cor- 
rections, G is related to K in a simple manner: 

G(R,pco) = — — 

K(R.pco) 
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Relative to Eq. (79) this appears to be almost the standard result if one were 
not in the mixed representation and if account had not been taken of the "ex- 
traneous coordinate" present in a non-uniform system. Work will always be 
done in the zero- order in the realization that this linearizing approximation 
can, in principle, be corrected by considering higher order terms in the 0 
operator. However, it is believed that 0 represents an asymptotic expansion 
which unphysically diverges in higher than zero- order terms. 


To extract physical quantities such as electron number densities, the 
following relations are needed. In analogy with Eq. (13), the mixed repre- 
sentation charge density is given by 


n(R,p) 


1 

2-n-i 



doo G(R, pco) 


(81) 


True physical charge density in coordinate space is given by 


n(R) = — ^-5 J d 3 p n(R,p) 
(2*r P<P F 


(82) 


so it is seen that n(R,p) is similar to a charge density in momentum space. 
If one uses a formal identity to write the Green's function as 

G(R,pco) = ■ ^ cj>(R) — , (83) 

u> - h(p) go - h(p) gj - h(p) - cj>(R) 

then one can see that the first term on the right-hand side gives the initial 
charge density when put in Eqs. (81) and (82) whereas the second term gives 
the screening charge density: 

6 n(R) = 5- f d w f d 3 p' — — 4 ><R) — 1 — 

2 Tri ( 2 ir ) 3 J p<p_ - h(p) oo -h(p) - 4 >( R) 

(84) 

= 7 f dcj / d 3 p / d 3 p' <|>(R) - 5 ( P^P' . - ^ 

2iTi ( 2 ir) J P<P F P ,:> P F w - h(p) w - h(p') - tj>( R) 


where q is the momentum transfer in the scattering event described below. 
The screening process may be envisioned in the following manner. An inci- 
dent electron of momentum p comes upon the impurity and scatters off of it 
transferring momentum q. The outgoing momentum is p’ 1 . The electron 
can scatter directly off the impurity, or through intermediate state many- 
body processes involving excitation of electron-hole pairs. Within the RPA, 
the lowest order impurity scattering processes are shown in Figure 10. The 
sum of all possible RPA processes is equivalent to a single scattering from 
a screened potential. Thus one can write: 
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Figure 1 0 - R epr e s entative diagrams included in static 
RPA screening of an impurity. 



1 iq. R 
e(q,o) e 


(85) 


Scatterings which involve more than a single momentum transfer to or from 
the source, such as those shown in Figure 11, are not considered, as is the 
customary procedure (ref. 57). If one approximates the outgoing propagator 
as a free-particle Green's function and puts Eq. (85) into Eq. (84), then 



Figure 1 1 .--Representative diagrams not included in 
static RPA screening of an impurity. 
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with the Pauli principle restrictions implied. Doing the p' and u> integrations 
results in 
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The p integration is done according to standard procedure (ref. 6) and yields 
an integrand which is simply related to € . Hence: 


3 

6n(R) = — ~ 

(2u y 



( 88 ) 


where now q and R are given in units of Fermi momentum. Equation (88) 
is further simplified to a one-dimensional integral by performing the straight- 
forward angular integrations. Thus: 

3 °o 

^=3 Unn-0 <89) 

cn q 

which is the result first obtained by Langer and Vosko (ref. 19). The charge 
density is given in units of 

3 

P F * 


Numerical examples appear in section VI, when bulk impurity screening is 
compared with surface impurity screening. 

The derivation of simple RPA screening has been covered for two 
reasons. The first was to see what really was hidden in the standard deriva- 
tions. One thing which is never cited is the assumption equivalent to using 
the zero order 9 operator. If one is interested, the formal machinery to 
compute higher order corrections exists. Another interesting observation is 
the equivalence of the pseudo-Hartree approximation to the non-local part of 
the mass operator and the assumptions implicit in a linear screening theory. 
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The second reason for going through a Green's function screening theory 
in which expansion is made in terms of screened rather than bare Coulomb 
potentials is that the present formalism is more general than the LV self- 
consistent field approach. A wider class of problems can be treated within 
the context of a single formalism developed here. In the next section the 
case of a surface impurity will be treated. 
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VI. IMPURITY SCREENING IN THE SURFACE REGION 
OF AN ELECTRON GAS 


Generalities 

Since the theoretical ideas- necessary for considering two different sorts 
of inhomogeneity problems, namely, surfaces and impurities, have been de- 
veloped thus far, it seems natural to attempt a coupling between the two, 
that is, impurities in the surface region. 

This problem should be of great interest for various reasons. First, it 
is of intrinsic interest in the general theory of inhomogeneous electron gases. 
A more practical reason is that this model system is precisely a model for 
ionic adsorption on metal surfaces. The role of ion screening in the surface 
region has only recently been gaining recognition (refs. 40-44), although a 
thorough theory has yet to be produced. This topic is further developed in 
this section. Previous adsorption theories have treated the metal or elec- 
tron gas in a very phenomenological sense; that is, it has been assumed that 
the ion-metal interaction would have basically the same functional form as 
a classical image force, even at short distances where field penetration and 
bulk screening occur (ref. 58). This section will show why these theories 
have been so successful in spite of their lack of rigor in treating the metal 
or electron gas. 


Theoretical Considerations 

The idea to be contemplated is this. A point impurity of charge Ze is 
placed on the surface of a semi- infinite electron gas. As was pointed out 
in Section IV (Dielectric Function), the response of the semi- infinite gas to 
an impurity of charge Ze placed at z = 0 is exactly the same as the response 
of the infinite sine wave electron gas to an impurity of charge 2Ze at z - 0. 
The dielectric function for the infinite sine wave gas is known, so such 
things as screening charge densities could be calculated. 

For various computational reasons, one is better off working in the co- 
ordinate representation initially since it is not yet clear how momentum 
conservation holds in this coupled- inhomogeneity problem. The starting 
point is the configuration space analog of the Dyson equation, Eq.(10), for 
the exact, interacting electron Green's function in the presence of an im- 
purity. 

G(r", J r';t) = G^rVr'jt) + J d 3 "r" G^rV^'jt) V (T") Gfr",''r';t) . ( 90 ) 

As mentioned in Section V, when computing charge densities, the first term 
on the right-hand side gives the unperturbed charge density whereas the 
second term gives the screening charge. If, as in the case of the bulk im- 
purity, one considers screening to occur within the context of a scattering 
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event, and if one considers processes involving only a single momentum 
transfer from the bare impurity, then G can be replaced by G Q within the 
integral of Eq. (90). The screening charge density is thus determined 
through knowledge of 


lim _^SG(r, r';t) 
"mr 1 
t—0 



G o (?,7"-,t) Y(r" ) G q (>,■?';!) . 


(91) 


The screened impurity potential is given by 



(92) 


where the factor of 2 comes from the image method and e g (q, 0) is the sine 
wave dielectric function, Eq. (66),whose z dependence is discussed later in 
this section. For completeness, the formula for € s (q, 0) when the integrals 
of Eq. (60) are carried out is cited: 


€ s (q, 0) - 1 
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In the above equation, a Q = the Bohr radius. 


The non-interacting Green's function is written inaneigenfunction ex- 
pansion as 
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Using the eigenfunctions given by Eq. (62), Eq. (94) becomes 
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If one writes the sine wave product in G 0 (r", r') as a sum of products of ex- 
ponentials and combines Eqs. (91), (92), and (95), the result is 
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Owing to the perfect symmetry or asymmetry of all the functions of k z and k^ 
in the integrand, it can be shown that appropriate variable transformations 
of the form k z — -k z in each of the four z' exponentials can be n ide. These 
simplify things greatly. When this straightforward procedure is followed, 
all four r' integrals take exactly the same form and give a result 
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which guarantees momentum conservation at the vertex, 
reduces to 


Thus, Eq. (96) 
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The product of sines is written as 
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-ez z -e z z I. 

In the physical system considered, an inversion symmetry about the z - 0 
plane exists. Thus the theory is invariant with respect to the transformation 
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z — ► - z, provided one keeps only the symmetric part of any oscillatory ex- 
ponential in the end results. Thus the sine function can be written in an 
equivalent form: 


r 

sin k z sink z 
z z 
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e - i( k 
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At this point, since an interest in charge densities given by Eq. (13) 
exists, the energy integral is performed. This step, together with Eq. 
(99) put into Eq. (97), gives: 
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Performing the k’ integration to use up the delta function yields: 
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Considerable insight into the problem can be gained by contemplating the 
structure of Eq. (101), which is an exact result so far. The first term is 
exactly the result of Langer and Vosko except for the fact that the non- 
spherically symmetric dielectric function appropriate to the sine wave 
electrons is being used. This term will present no computational problem. 
The second term has an extra oscillatory term in it as a result of pure 
surface effects on the incident and scattered electrons. One might hope 
that some physical simplifications will allow evaluation of this expression 
in the interesting asymptotic regions of z small and large. Note that the 
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second integral has virtually no hope of being exactly evaluated since the 
standard methods of treating the Pauli principle restrictions can be meaning- 
fully invoked only for spherically symmetric k integrals. However, one can 
still proceed with his exact treatment, even if only in a formal sense. The 
idea is to remain exact for as long as possible before contaminating the 
theory with approximations needed to obtain numbers. 

There exists a form of the mean value theorem which now is of use 
(ref. 59). Theorem: If there is an integral of the form 

b 

I = J f(x) g(x) dx, 
a 

with f(x) and g(x) well behaved, then it can be proved that there exists a £ 
with a S £ < b , such that 

b Id 

J f(x) g(x) dx = f (i) J g(x) dx. 

a a 

— Z ilc z 

For this particular case, one takes f(x) = e z Thus Eq. (101) 

can be written 


6 n^) 




y(q) 

€ s (q: °) 



The k integration, as in Eq. (87), is simply related to the spherically 
symmetric RPA dielectric function by v(q) I (q, k) = e (q, 0) - 1 . Performing 

this integral yields: 


5n (?) = f d 3 q-4^r e iq * r (e (q) - 1) (1 - e _2l ^ Z ), (103) 

(2 tt) J s [q> 

which is still an exact result. 

Now the more difficult aspects of the theory must be faced. Two prob- 
lems are present. First, there is no simple way of determining £ and, 
second, £ depends parametrically on both q and z. Fortunately the physics 
of the situation can help some. In reference to Eq. (101), it was mentioned 
that the oscillatory term of the second integral was a result of the surface 
boundary conditions. This effect is washed out within the order of a 
Fermi wavelength. Mathematically this can be demonstrated easily. For 
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large z, the exponential oscillates rapidly in momentum space, and thus the 

integral on k — 0 by virtue of the Riemann-Lebesque lemma* 
z 

Lemma: 


If J f(k) dk <oo 3 

then 


lim 


z— la 


rge J 


f(k) e lzk dk— 0. 


Thus the exact long-range screening charge density is given by 




d 3 q-L= (€(q) - 1) e 1 '!' r . 




(104) 


It is fortunate that the long-range part can be obtained in an asymptotically 
exact form for this enables us to determine the effect of the surface on 
impurity Friedel oscillations with no approximations. 

The short-range behavior is another problem. In any approximation 
scheme, the most important condition to be maintained is that at the sur- 
face z = o, the screening charge density is zero, for in a linear dielectric 
response model there certainly can be no induced charge where originally 
there was no charge. Again, in reference to Eq. (101), one might expect 
that in the near region things would not be too dependent on the accuracy of 
that is, for z small, the exponential varies slowly relative to other 
factors in the integrand. Thus an average value of k within the integration 
range 

(0< k z < P F ), 

say, 



is arbitrarily assigned to £. This is a good approximation as z— 0 and 
hopefully is an adequate one for intermediate ranges where z = 0(p^,). 

With this, the first approximation, and the replacement of the exponent by 
the symmetric cosine as discussed previously, Eq. (103) in the near region, 
takes the form: 
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d 3 q-±z (€ (q) - 1 ) e iq ' r 
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(104) 


or with an elementary trigonometric identity: 

(2 ' \l2 


6 n (r) = 
s . r 


"■’(W 


2 . 2 f**F" A I ,3 1 , , . .. iq-r 

d ( 6 (q) “ 1) e H . (105) 


Equation (105) is easily interpreted. The first thing to notice is that the 
surface boundary condition at z = 0 is satisfied. Second, one can see 
through the structure of this result. The terms within the integral are just 
the same terms as in the long-range part (Eq. (104)). These terms reflect 
the many- body impurity Friedel oscillations modified by the surface dielec- 
tric function. However, in the surface region, z = small, surface boundary 
condition density oscillations must play a role. Hence, the presence of the 
multiplicative "surface sine" term is expected. The resulting charge 
density configuration has the form of sine wave screened impurity Friedel 
oscillations superposed on surface density oscillations, a not unreasonable 
result and incidentally the sort of result Kohn and Sham (ref. 30) contemplated 
in problems of impurity screening in inhomogeneous electron gases. 


One word of explanation as to the nature and motivation for the present 
method of approximation is in order. Although it may be true that the 
asymptotic results could have been surmised without going through the for- 
mal development since Eq. (101), it is felt that the formalism has been 
useful for various reasons. Equation (103) is an exact result. Although the 
author has not been able to conceive of a truly convincing general approxi- 
mation for the quantity £, he feels that there probably is a reasonable pro- 
cedure such as a modified stationary phase or steepest descent method by 
which £ could be obtained. If this procedure could be discovered, then the 
present theory up to Eq. (103), which is the end point of the theory, would 
be immediately applicable with no modifications. Another possible way of 
determining an approximate value of £ would be to convince one's self that £ 
is not very q- dependent, but is only a function of P-p* Then one could take 
the volume integral of 6 n(r), set it equal to 1 , and see what value of £ is 
needed to guarantee self-consistency. This is left as an excercise for the 
future . 


The non- isotropic surface dielectric function will now be treated. From 
Eq. ( 66 ), it is seen that € (q) can be written as 


« s (q) = 2 <(q) + 1 e(q) 6 q ,0 


(106) 
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However, it must be remembered that this is the dielectric function at 
z = 0 only. In effect, the above should really be writteii as 

6 s (q) = 6 s (qi z = o). 


As z increases, the dielectric function must approach the volume limit, 


lim z-large e S W z) = ' (q) ' 


Folowing lines similar to Stern (ref. 42), one calls upon an ansatz for 
dealing with the z dependence by allowing the surface effects to decay ex- 
ponentially into the metal while the volume effects increase in the same 
manner. With the ansatz, Eq. (106) becomes: 


€ g (q;z) = \ € (q) + \ * (q) (1 - e" qZ ) + ~ e (q) e" qZ 6 (107) 

4 z , 


which upon rearrangement is 


,(qjz) - € (q) + j e (q) e 


-qz 


1 ). 


(108) 


Equation (108) clearly satisfies both the small- and large-z-limiting con- 
ditions. 


The inverse dielectric function when taken inside an integral on q, that 
of Eq. (103), is of concern. Using the operator identity 


1 = 1 _ J_ b 1 

a + b a a a + b 
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for e (q;z) in Eq. (103) with a = e (q) and b = e (q)e 
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which can formally be written as 



1 

*(q) 




q 2 t o. 


(UO) 


Equation (110) is a valid representation of Eq. (109), only when e (q;z) 
appears within an integral, as it does in the present case. With arguments 
as given in section IV, it is seen that the q ^ 0 restriction on the second 
term provides no difficulties since the integrand is well behaved for q^ = 0- 
Consequently, one obtains the simple form for the inverse dielectric 
function: 


1 _ _J_ / 2 

€ g (q;z) “ € (q) \ 2 _ e -q z 

With Eq. (Ill), the expression for the surface impurity screening charge 
density (Eq. (103)) becomes: 


(HI) 


6n (r) 




lq ♦ r 


( 1 


( 112 ) 


If all length and momentum, variables are now written in units of Fermi 
momentum, and the angular integrals are performed, then the asymptotic 
limits of the screening charge density become: 



o 


These terms represent the Langer and Vosko result with two modifications. 
The multiplicative sine factor merely says that surface density oscillations 
must be superimposed on the impurity quantum density oscillations. The 
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z dependent term within the integrand is also a departure from Langer and 
Vosko's result. This is the term which couples what happens to the 
spherically symmetric response as a result of the constraints in the z 
direction. 

In the next section, the numerical evaluation of these integrals is 
presented, the results discussed, and experimental suggestions made. 

Numerical Results and Discussion 

Equations (113) and (114) have been evaluated numerically for several 

different electron gas densities along various radial lines emanating from 

the impurity. The results are displayed graphically in Figures 12 through 

15 for values of r = 1.5, 2, 3, and 4. 
s 

Much can be said about these results of screening charge as a function 
of vector position from the impurity. One should first observe that these 
results are quite understandable. At R = 0, for any Q , the initial charge 
density is zero because of the infinite barrier. Thus, in the linearized 
model, there will be no polarization charge where initially there was no 
charge. As the weakly screened field penetrates into the gas, polarization 
charge builds up as the unperturbed electron density increases. Further 
into the electron gas past the point where the maximum screening charge 
density is seen, the screening becomes more effective and the polarization 
charge density decreases. Long-range Friedel oscillations, plotted on an 
expanded scale in the insets, are also present, although their structure is 
more complicated than volume impurity oscillations. This results from 
the fact that a phenomenon involving interferences between two diffraction 
patterns of different symmetries, that of the surface and that of the im- 
purity, is being considered. 

For a quantitative comparison with volume effects, the results obtained 
by Langer and Vosko for an impurity in a homogeneous electron gas of 
density r = 3 are also plotted in Figure 14. It can be seen that the surface 
impurity screening is of a much longer range than the volume impurity 
screening, as it must be. In the surface problem, the low-density region 
cannot provide enough electronic charge to do much screening. Hence the 
screening charge is less than in the volume case. To make up for this 
deficiency, as the relatively unscreened surface impurity field penetrates 
the metal, it must induce a larger long-range polarization charge than in 
the volume case, so that the impurity may be totally screened at large 
distances. Consequently the surface impurity charge density curves cross 
over the volume curves. 

As was mentioned, the structure of the oscillatory curves is more 
complicated for the surface impurity, at least in the intermediate region 
4 <Rpp £,10, where chaotic things are happening. As this chaos settles, 
the longer range Friedel oscillations, say for Rpp £,10, become well be- 
haved in the sense of volume oscillations. Still there are significantly 
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Figure 15 - Surface impurity screening charge density, 


r =4.0. 
s 
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longer range effects for the surface impurity than the volume impurity as 
suggested in the preliminary theory of surface impurities (ref. 43). This 
observation could be of experimental interest for a number of reasons. 

First, standard magnetic resonance and Knight shift experiments are 
measuring effects occurring within the skin depth of metals. Since the con- 
duction band electron polarization around ion cores or impurities in the 
outer skin is different than in the bulk, measurable events, such as impurity- 
impurity interactions, may be different than those inferred from standard 
theories using the customary long-range polarization charge density 

Sn ~ (cos 2p^,R / R^. 

Second, the existence of polarization properties different on the surface 
than in the bulk should offer encouragement for devising new resonance- 
type surface experiments to replace the rather dreary assemblage of pre- 
sently existing experiments. It should also be noted that in the case of an 
adsorbed partial monolayer, a distribution of surface impurities, there may 
be significant impurity- impurity interactions caused by the intermediary of 
the polarizable electron gas. Put another way, the long-range surface 
Friedel oscillations, for large angles, of one impurity may interact with the 
oscillations from a nearby impurity, thus giving rise to an effective inter- 
action which could not be inferred without considering the many-body 
screening effects. It can be seen that the problem is much more subtle than 
would be guessed from the various oversimplified theories which exist. 

With regard to the numerical results, a few more comments are in order. 
As anticipated, it can be seen that the screening charge increases more 
slowly as the transverse component of the radial line increases, since the 
unperturbed electron density doing the screening increases more slowly. 

The following structure can also be observed in all the curves. A 
maximum in the screening charge density is seen at Rp F ~ 1. 1 and a second 
rise in charge density at Rp 5-6 corresponding to the first Friedel 
oscillation. However, the amount of screening charge in the first lobe 
decreases as the electron gas density decreases. For instance, if r =1.5, 
then 


An(Rpp= 1.1) = 0. 0425 for Q - 0 
whereas if 

r g = 4, An (Rp F = 1. 1) = 0. 025. 

This means that in the lower density gas, more of the screening, and thus 
induced charge, appears further into the metal. For this reason, the 
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classical image charge is not applicable for honestly describing surface 

impurity effects because the screening charge exists in the metal and not 

on the surface. This, of course, is a manifestation of the fact that the 

electron gas or real metal is not a perfect conductor as must be assumed 

if the concept of an image charge is to contain any real meaning. In the 

case of a perfect conductor, r =0. The total screening charge will appear 

within the first lobe with a maximum at Rp_ ~ 1. 1. But since 

F 

P F = , lim _ R~ r — 0 . 

* r r -~0 s 

s s 

Hence the screening charge appears either at the impurity site with R = 0 
or in the 6 - 90- deg curve, in either case at z = 0. This is then a true 
surface polarization charge which is a consequence of the necessary 
conditions for having a perfect conductor. As the density is decreased to 
real metal densities, then screening effects and their consequences must be 
considered from the volume polarization point of view. The fortuitous 
success people have had using the image force and surface polarization 
concepts in absolutely no way says that the approach is correct. 

The present results may be applied to existing theories of ionic 
adsorption in the following manner. Essentially one is interested in the 
dipole moment formed by the point ion core or impurity and the resulting 
polarization charge since the change in electron work function is given by 

A d> - 2v cr M 
' e 

where cr is the density of dipoles and M is the dipole strength. As stated in 
Eq. 23-1/2 of the author's atom -metal interaction theory (ref. 8), the dipole 
strength is M = q (s + e (s) ) where q = the ion charge, s = the ion- 
surface separation which could be taken to be the ionic radius, and "e(s) is 
some effective location for the polarization charge of the positive ion core 
and has a numerical value greater than s for small s." The assignment of 
the value of the ionic radius for s is dubious at best. However, since this 
seems to be a convenient way of assigning parameters to the size of the 
adsorbate, the idea is introduced. 

A value for € (s) can be obtained from the calculated screening charge 
densities as follows. Approximate calculations of the dipole component 
of the polarization charge distribution 

M = q < z> = / z 6n (R) R 2 dR dQ (115) 

have shown that a reasonable representation of the calculation is given by 

<z> ~ — * — = .96r A. (116) 

P F s 
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Thus it is apparent that the strength of the effective dipole moment of an ion 
on a real metal surface is quite dependent upon the substrate properties, the 
free electron density. This observation has not previously been subjected to 
quantitative analysis, and, with only a few exceptions, has not even been 
qualitatively recognized. To apply these ideas in real calculations, one must 
obtain the density of s- and p-like conduction band electrons, those electrons 
which participate in the screening. If one is dealing with noble or transition 
metals in which there is great hybridization between the free electron- like 
s and p bands and the tight- binding d bands, one must not include the electron 
densities coming from the d bands since the d electrons do not contribute 
much to the true screening charge density. 

In a typical system of experimental interest, such as cesium on 
tungsten, sample parameters for theoretical calculations might be 

o 

s = r . = 1 . 3 A 

i 


and 


2<r s <3. 


With this choice of fundamental atomic and metallic parameters, the dipole 
length would be: 


o 

— — = (r. + . 96 r ) — 4 A. 
q l s 


(117) 


This is in accord with work function data. The point to be made here is that, 
in the calculations of o reference 58, it had to be assumed that the dipole 
length was at least 4 A or the atom-metal separation at least 2 A in order 
that meaningful results might be obtained. This separation, however, is 
quite a bit greater than 1. 3 A which must be used in an image force picture. 
Now it can be seen why the 4- A dipole length assumed is the correct value 
and not the 2. 5-A value required in a classical image force model. The 
general rule is this: If the ion-metal separation becomes as small as a 

Fermi wavelength, then the average position of the screening charge will be 
deeper into the metal than the ion-metal separation. In this case, perfect 
conductor image forces have little quantitative significance. With the 
simple approximate relation given by Eq. (116), it is imagined that more 
realistic data correlation with relevant properties of both the ion and metal 
systems will now be possible. 

One final point to be made is concerned with the general trend of work 
function changes with s and p electron densities. From Eq. (116), it would 
appear that the work function change would be greatest for an ion adsorbed 
on a metal with the lowest electron density (r g largest). However, another 
competing trend must be considered. In general, the bare work function of 
a metal decreases with decreasing electron density. Thus, the position of 
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the shifted and broadened atomic level will overlap more of the occupied por- 
tion of the conduction band. Consequently, the atom-metal bond will become 
more metallic and thus less ionic in character (ref. 58). In this case, the 
effective charge on the adsorbate is reduced and the actual dipole moment 
M = qeff( r i + <2 ^) will decrease. 

There are thus two competing trends. As the electron density is de- 
creased, <z> increases but q e ff decreases. Hence there exists an optimal 
electron gas density for the greatest work function change using a particular 
adsorbate to specify the value of r^. Put another way, it may be that cesium 
would produce the greatest work function change in one metal wheras potas- 
sium would produce the greatest change in some other metal, solely as a re- 
sult of the different electron gas densities in the two metals. 

It is hoped that considerations along the lines outlined in this section 
will be taken by workers doing numerical calculations of electronic surface 
properties for it seems that if such considerations are not given, then there 
can be little correspondence between what is calculated and what is occur- 
ring in physical reality. 
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VII. INFLUENCE OF MANY- BODY EFFECTS 
ON THE SURFACE POTENTIAL 


Generalitie s 

A perhaps mysterious claim was mentioned in the introductory section 
in which an allusion to the necessary connection between the surface poten- 
tial barrier and many-body effects was made. A demonstration of this as- 
sertion and how existing theoretical treatments of the surface barrier can 
be improved upon by using the theoretical ideas set forth in the preceding 
sections will be presented in this section. 

Basically, the idea of a surface barrier infers that deep within a metal 
an electron propagates through a partially non-local potential generated by 
the many-body medium of the electron gas, plus the periodic pseudo-poten- 
tial of the lattice. As shown in section III, the periodic lattice does not 
drastically alter the many-body effects. As the electron travels about, it 
undergoes ene rgy- reducing processes of the exchange and self-energy type. 
For the most part, these terms do not depend upon position because of the 
translational invariance of the system. However, thinking in Thomas-Fermi 
terms, one would expect that if the electron gas density varied, then the ex- 
change and correlation effects, depending directly upon the electron gas den- 
sity, would also vary. In the simplest possible illustrative example, con- 
sider the average Hartree Fock free electron exchange energy = -0.916/ r g . 
As the density decreases, r s increases and thus the reduction in total energy 
caused by exchange decreases. Hohenberg and Kohn (ref. 29) have given a 
general theorem which proves that the energy and density are uniquely re- 
lated. Although it is the author's opinion that the given proof of this theorem 
does not apply to the density variations in the surface region, „ he believes 
that the physical idea is basically sound. As the propagating electron comes 
into the surface region, the electron density of the field within which it is 
propagating decreases. Thus one would expect the energy lowering exchange 
and self-energy effects to be less pronounced. Consequently, the electron 
"potential energy" (a loose definition of that part of the total energy not in- 
cluding kinetic energies) increases in the surface region. However, the 
electrons are in energy eigenstates so the kinetic energy must go negative, 
corresponding to the classically forbidden region. Therefore binding occurs 
as a result of the many-body effective potential barrier. An interesting 
item to contemplate in this line of thinking is the following. Since the bar- 
rier is caused by many-body effects, such as self-energy processes, the 
barrier height depends on the momentum of the particular electron incident 
upon it; that is, an electron on the Fermi surface which is normally incident 
upon the metal surface sees the outside vacuum potential at a height 4> e# the 
electron work function. However, an electron with kinetic energy p^/ 2m - e 
normally incident does not see a height 4> e + e because its momentum is dif- 
ferent. The momentum- dependent energy reducing many-body effects it in- 
duces are different. Consequently, when referring to such phenomenological 
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constructs as electron work functions, one should use a relation of the form 
^ e “ ^ed^) with the condition 4 > e (pp i z ) = ^p. s ua i where the true barrier, 
minus the Fermi energy which a normally incident Fermi surface electron 
sees, has been arbitrarily set equal to the phenomenological barrier minus 
the Fermi energy. This sort of phenomenon might have quantitative signifi- 
cance if one derived the correct Richardson- type equation for thermionic 
emission. In the derivation it is necessary to do phase space integrals of 
the form 


/ d 3 p e- *' kT f(p. T, 
which should be of the form 


kT f( - T) . 

From the relative success of the standard Richardson equation in interpret- 
ing experimental data, it is probable that these corrections are not of first 
order importance. However, as the quality of the experimental situation 
improves, these observations may be of more than academic interest. 

There are other far-reaching motivations for studying the nature of the 
surface barrier. Periodic deviations in the Schottky effect have been the- 
oretically studied by Cutler and Gibbons (ref. 61). They demonstrate the 
strong dependence of these deviations on the reflection coefficient for a 
metal electron incident upon the surface barrier. It is well known from ele- 
mentary quantum mechanics that the reflection coefficient is extremely sen- 
sitive to the form of the barrier and its higher order derivatives. D'Haenens 
and Coomes (ref. 62) have displayed the problems one gets into by oversim- 
plifying the nature of the surface barrier. In attempting to explain their ex- 
perimental data, they use a Sommerfeld box model for the metal with no re- 
gard to many- body effects. Using a perfectly square barrier to represent 
the surface, they require bandwidths which are roughly twice as large as in 
reality in order to achieve agreement with experiment. This of course is 
an entirely unacceptable situation and the idea of the square barrier must 
be discarded in calculations requiring reflection coefficients. An adequate 
replacement is required. 

Still another field of research should be concerned with the nature of 
the surface potential barrier, the field of photoemission studies. Spicer 
and co-workers (ref. 63) have developed experimental techniques for doing 
photo studies to a high degree of accuracy. From measurements of energy 
distributions of emitted electrons, they infer bulk density of states and at- 
tempt to explain certain optical properties of solids. Decimal point accu- 
racy assignments are given to particular interband transitions, based on a 
simple theoretical model for the photoemission process. It seems that a 
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more complete study of the surface barrier and its effect upon an incident 
photo -excited electron would help to understand further the work of Spicer. 

This problem was considered initially by Bardeen (ref. 33), Juretschke 
(ref. 34) performed a calculation of the unscreened exchange potential as a 
function of the distance from the surface which demonstrated both the exis- 
tence of the barrier and the fact that the potential is not monotonic but has 
oscillatory structure. Loucks and Cutler (ref. 35) refined this calculation 
by using an exchange potential shielded by a spherically symmetric screen- 
ing function and including rudimentary correlation effects within the context 
of the plasma oscillation theory. Here the Loucks and Cutler theory is ex- 
tended by considering collective effects more appropriate for an electron 
gas with a surface, although still in a fairly rudimentary manner. However, 
a formalism is set up which would enable one to do a more exact calculation 
if this is felt necessary. As mentioned earlier, because of the present state 
of the art on the experimental side of the fence, such undue effort does not 
seem justified. 


Theoretical Considerations 

The objective of this section was to calculate the self- energy for a prop- 
agating electron in an electron gas with a surface. The difference between 
the self- energy in the surface region and that deep within the interior gives 
a contribution to the surface potential barrier. 

The physical picture can be stated in the following manner. Initially 
one starts with a propagating plane wave electron. If it happens to enter the 
surface region, the density fluctuations it induces will be different from those 
in the interior. Consequently, in the case of electrons with kinetic energies 
less than the Fermi energy plus phenomenological work function, scattering 
occurs back into the metal. The nature of the surface is seen through this 
picture: an incident plane wave electron scatters off of a density fluctuation 

it induces in a much different manner than in the bulk. The scattering is 
such that an incident electron of momentum k z will have outgoing momentum 
of -k z as a result of scattering from the surface potential, V g (z). Thus the 
sine wave functions previously used in surface calculations are really scat- 
tering states, the sum of an incident plus scattered wave. An attempt to 
calculate the form of V s (z), the effective potential induced by an incident 
plane wave, will now be made. 

The method utilized in the calculation is now outlined. As discussed in 
section II, the self- energy of an electron propagating in the many- body field 
can be written as 


S (p, p°) 



G(p-q) r (p,q) D(q) 


(118) 
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with 


D(q)=^L 

€ (q,q ) 


An equivalent statement of Eq. (118) is: 



(119) 



within which one could introduce the density fluctuation propagator: 



( 120 ) 


In the first integral of Eq. (119), if one sets G = G 0 and r _ 1, then this is 
just the unscreened Hartree-Fock exchange contribution to the self- energy. 
The second term in Eq. (119) involving the propagation of the density fluc- 
tuation can also be split up into two parts -- one involving density fluctua- 
tions arising from excitation of virtual electron- hole pairs and the other in- 
volving the excitation of a virtual plasmon. The procedure for effecting this 
division is discussed by DuBois (ref. 6). Here his end result is simply 
quoted. It says that 


The first term, which takes into account pair excitation, gives rise to modi- 
fied short-range interaction effects. The point of view is adopted that if term 
number one was dealt with in a rigorous manner, then one could derive the 
screened exchange self- energy. Then, if one devised a method of modifying 
the bulk theory to a surface theory, the end result of a proper treatment of 
the Hartree-Fock plus screening terms would give the screened exchange 
potential calculated by Loucks and Cutler (ref. 35). 

Within the context of this point of view, the only thing left to consider in 
the self- energy is the second term of Eq. (121) which is related to the plasmon 



2[E o (k + q) - E Q (k)] 


78 


part of a density fluctuation. In this term, the 6 step function is introduced 
to indicate that the plasmon is a well-defined excitation only for momenta 
less than some critical momentum q c where the relatively flat plasma dis- 
persion curve merges with the pair continuum shown in Figure 16. 



Figure 16. Elementary excitation spectrum in an electron gas. 

It will be convenient to introduce a parameter (3= q c /pp. 

It might be expected that since the q = 0 collective oscillations in the 
surface region have a frequency to g = Up/ as compared to w in the in- 
terior, a difference in plasma self-energy effects between the volume and 
surface region could be categorized by properly taking this effect into ac- 
count. The problem is that in the surface region, the surface plasmon has 
its momentum parallel to the surface only. Thus, in the surface region, 
limq^ — -o w s (q”) = w p/ , whereas the collective oscillations with momen- 
tum normal to the surface go as lim^ — q ^ (q) = Up, To perform the surface 
calculation with absolute correctness involves using a non- spherically sym- 
metric plasmon dispersion relation to account for the fact that surface plas- 
mons propagate parallel to the surface only. In the present approximate 
treatment, it may be assumed that the dispersion relation is symmetric in 
the sense that an electron in the surface region excites a virtual plasmon of 
frequency cjpA/2 independent of angle, whereas an electron in the bulk ex- 
cites one at w p- This is an improvement upon the Loucks and Cutler treat- 
ment in which the relevant plasmon contribution was taken at a frequency of 
oj p independent of position. 
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With these assumptions, one can now proceed with the mathematical as- 
pects of the theory. The plasmon part of the self-energy of Eq. (119) can be 
written as 


2 (P, P°) 
P 




Dp(q) G(p-q) 


if, by use of Eq. (121) in (119), one identifies: 

2co-*- 

D ( q ) = — 2 2 . 

^ co - OW+ 1 -n 

q ' 

which is a standard Boson propagator for the plasmon where 


( 122 ) 


(123) 


v('q) 

8c (q^co^)/ 3u 


(124) 


is an effective coupling between the propagating electron and virtual plasmon. 
The mathematical aspects of this problem, a propagating Fermion emitting 
a virtual Boson in a self-energy process, are identical with many other prob- 
lems in physics such as the polaron problem (ref. 64), He^ coupled to phonon 
modes (ref. 65), and hole coupling to plasmons (ref. 10). 


In the long wavelength limit, the dielectric function is given by 




1 - 


2 

co 


s 



Furthermore, it is a reasonably good approximation to take the surface plas- 
mon as being dispersionless; that is, w q~= <^> s = constant. With these assump- 
tions, Eq. (124) becomes: 



Equations (123) and (124) in Eq. (122) give: 


(125) 


/ 4 

“^4^2 ~ 2 \ G( P~ C 1 ) (126) 

(2-rr) "q" to - co^ + i T] 
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The next item on the agenda is to consider the appropriate Green's 
function for the incident electron G(p-q). This can be formally written as 


G(p-q) - (127) 

+ i Ti 
p-q ' 

with the following explanation. If one takes Ep_q to be the energy of an elec- 
tron deep in the metal, and if one then sets the exchange potential at z = 0 
as the zero, then 

E^- — = ( P~ T g; ) + { V . - V (0)1 

p-q 2 m | int x J 

results. In this expression V^ n £ is the total internal potential which includes 
both screened exchange and correlation effects. On the other hand, V x (0) is 
only the screened exchange part of the surface potential since it is precisely 
the correlation energy part which is to be calculated. Note that V^ n f. and V x 
are negative numbers. For brevity. 


Z 

int 


V 


int 


- V x< 0) 


may be written. 

Some further approximations may now be introduced. These approxima- 
tions will limit the accuracy of the calculations to within a factor of 2 or so. 
But this is sufficiently accurate to learn whether a more detailed calculation 
is worthwhile. The approximations are this. First, one can assume that 
q << p in E°_q. Since one is only considering q < q^ for plasmons, the 
first assumption does not completely tax the credibility gap. Next assume 
that the dominant processes occur for electrons near the Fermi surface and 
thus replace p ~ Pp. With these assumptions, the electron Green's function 
becomes 


G(p-q) a 


128) 


2m 


Z + i r| 
int 


where Z is also taken as a constant independent of momentum, 
int 
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Equation (126) can now be written as 


"surf 


. „ 2 2 

1 4 ire co 

s 


if 

q<q„ 


d q dco _1_ 


(2 tt ) 


2 (go + co ) (co - co ) 
q s s 


co - 


2 

P F 

2m 


S int 


(129) 


The energy integral is performed around the contour shown in Figure 17 
which encloses the poles at 


2 

P F 

co = +co , ~ — + 2 . , 
s 2m int 


but not co = -co g . Equation (129) then reduces to 


co 


2 , = - 2 -it e 

surf 


/ 


d q 1 

3 2 

(2irT q 


co + 2. + j , 

s mt 2m 

q< q c 


(130) 


which is simply 


co 


J surf 


— k |5 
IT F r s 


2 

P p 

co + 2. A ^ 

s mt 2m 


(131) 


For the case of the dispersionless plasmon (3 , the surface plasmon cutoff 
is related to |3 the volume cutoff by 


fi - 2- 1 / 4 

r s 


J3 = . 84(3 


If a similar derivation were performed for the plasmon contribution 
to the self-energy deep within the metal, one would obtain: 


2 

2 i = k „ p, 

plasma it F r 


(132) 


the standard result. As an aside, it is encouraging to see that the field 
theoretic approach, coupled with the approximations that have been made, 
is equivalent to the Bohm- Pines theory. In the next section it is discussed 


82 



I 



how these results are used in calculating surface potentials, and numerical 
results, which are the Loucks and Cutler calculation with the surface plasmon 
correction, are presented. The relevance of this correction to the experi- 
mental situation is also covered. 

Numerical Results and Discussion 

As has been previously stressed, the theory presented in this section 
has been directed towards obtaining corrections to the Loucks and Cutler 
results by considering modifications to bulk correlation energy processes 
appropriate in the surface region. Mathematically these modifications came 
in two different ways. First, it was assumed that a propagating electron 
couples to virtual surface plasmons of frequency u / j'Z rather than co in an 
isotropic manner. Although this is not rigorous, ^the end result giveS in Eq. 
(129) would be rather insensitive to a fractional change such as 1/ ^/Z in the 
resonant frequency. Second, a new Green's function, Eq. (129), was used 
to describe an electron coming into the surface region from within the metal. 

At this point one may wonder how the z dependence of the self- energy 
was handled, since it does not appear in Eq. (131), nor in the derivation 
leading to it. The explanation is this. What has just been calculated is the 
self- energy appropriate to the z = 0 plane only. A calculation for a general 
point z = 0 could be done by using the image technique, but the labor in- 
volved would be forbidding. However, since interest lies essentially in the 
form of the effective potential as a function of z, the following approximation 
scheme is adopted. First, it is assumed that the oscillatory structure of a 
fully calculated E(z) would parallel the structure of the Loucks and Cutler 
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potential energy curves. Within this approximation, one can incorporate 
the present results into their calculated curves essentially by a rescaling 
of the energy in the following manner. In Figure 18 their results are re- 
produced. These curves represent one- electron potentials which include 
screened exchange plus a z- independent correlation energy. The z- 
dependence arises only from variations of the screened exchange. 

If one considers the theory developed in this section, then the following 
modification to Loucks' and Cutler's work is required. At z = 0, the 
plasmon self- energy part of the "potential" differs from the interior value 
by a factor 


0. 84 co 

s 

2 

P F 

co + 2 . + .T — 

s mt 2m 

In other words for any (3, at z = 0 the Loucks and Cutler result should have 
added to it: 


AS = 2 , - 2. , 

surf mt 


+ e 2 'kpp 


. 16 co - 2. , 

s mt 


2 

P F 

2m 


2 
p F 

co + 2. + ~ 

l_ s mt 2m 


(133) 


Thus, the procedure followed here is to take one of their particular curves. 
The interior potential determines the V = V( °o ) ordinate. At z = 0, re- 
scaling is performed in the sense that V(0) = V^ffl) + A 2 is required. 

Then with this new scale, values of V(z) for z arbitrary can be read 
from their result. These steps have been performed and the results drawn 
in Figure 19 with electron densities taken at the value for Na. A number of 
interesting observations can be made when comparing the original Loucks- 
Cutler curves of Figure 18 with the surface- corrected ones. First, it 
should be noted that for all practical purposes, both sets of curves with 
|S < 0.6 coincide. This is reasonable though for two reasons. Plasmon 
contributions to total energies are being considered. The lower (3 is, the 
fewer plasma modes there are. Thus one would expect low- p systems 
to have the smallest plasmon self-energies and consequently the smallest 
discrepancy between the original Loucks- Cutler treatment and the present 
one. Another measure of the importance of the plasmon self- energy term 
is how it stands in proportion to other energy- reducing terms, in particular, 
the screened exchange term. As Loucks and Cutler show, the magnitude of 
the screened exchange decreases with increasing (3. Thus as (3 gets larger, 
the plasmon self-energy process becomes as important as the screened ex- 
change. As a result, one would expect corrections to the surface plasmon 
self- energy diagrams to become more important. Hence, two distinctly 
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different reasons exist for anticipating the surface plasmon correlation 
energies, as calculated in this section, to be more pronounced in the 
domain of higher p's. Comparing Figures 18 and 19, one can see that 
what was expected did actually take place. The height of the repulsive 
barrier at z = 0 for p = 0.8 and 1.0 is significantly greater in the present 
theory than in that of Loucks and Cutler. But this is reasonable for 
several reasons also. In the first place, the correlation energy, which is 
always attractive at reasonable electron densities, becomes less attractive 
as the electron density is reduced. Thus one would expect an appropriate 
correlation energy contribution in the surface region, the region of re- 
duced electron density, to be smaller than in the interior. The present 
calculations bear this out. On the other hand, in the original Loucks- 
Cutler calculations, this type of correction is applied to the screened ex- 
change energy but not to the correlation energy. Thus, as p increases, 
their surface potential curves become quite flat, not becoming increasingly 
repulsive fast enough as the surface region is entered from the interior. 
Since the value of p for Na seems to be ~ 0. 8, the surface correlation effects 
are pronounced and consideration must be given to those effects before a 
theory of real world physical systems hopes to be reasonably complete. 

One final observation is concerned with the possible reflection co- 
efficients at these potential barriers. Since the original Loucks- Cutler 
curves had a relatively small slope, the reflection coefficient would be 
smaller than in the present curves. This change in the gradient might be 
of some importance when one constructs simple, analytical approximate 
surface potential functions based on the results shown in Figures 18 and 
19 {ref. 66). 
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VHI. CONCLUSIONS 


Thus far this document has shown a fairly quantitative understand- 
ing of selected examples of inhomogenous interacting electron gases. In 
fact, most of the problems that have been undertaken have been exactly 
soluble, at least in the asymptotic sense and within the RPA or some more 
appropriate plane wave many-body theory approximation. This would seem 
to be just opposite what would be expected when one looks at the difficulties 
the Kohn school has had in obtaining practical solutions to inhomogenous 
electron gas problems. 

As might be expected, there is, of course, good reason why the type 
of problems dealt with in this report have been amenable to rather neat 
solutions in spite of their initially ominous appearance. The fundamental 
problem with non- translationally invariant systems has been that, in 
general, nothing could be said about the meaningful conserved quantity 
usually associated with momentum in an invariant system. Thus, when 
Kohn starts theorizing on a slowly varying electron gas, the closest he can 
come to anything resembling states which are superpositions of momentum 
eigenstates, are WKB wave functions. The inhomogeneous nature of the 
electron gas is thus reflected in the different single particle WKB energy 
eigenfunctions. However, these functions are difficult in momentum space 
so he cannot Fourier transform his configuration space problem into 
momentum space and take advantage of any conservation laws which 
greatly reduce the computational complexity of the problem. But all is 
not bleak because the problems which he deals with are those in which the 
change in total electron density relative to the average density is very 
small over a distance of the order of a Fermi wavelength. Consequently 
he can avail himself of very truncated Taylor or density- gradient ex- 
pansions to take care of the spatial variations. As the density change be- 
comes more rapid, higher order terms in the expansions must be main- 
tained in order that there remain any correspondence between calculated 
results and physical reality. When this situation prevails, the Kohn 
approach, with density as a basic variable, becomes of decreasing utility 
for obtaining results which will be evaluated numerically and compared with 
experimental situations. 

But, as if there were some overall justice in nature, those examples 
of inhomogeneous electron gases which do vary significantly in the dis- 
tance of a Fermi wavelength have other characteristics which make them 
calculable also. A little reflection by the reader will show that most of 
the physically relevant, rapidly varying gases come about as a result of a 
localized boundary condition. Take for instance the cases investigated in 
this report. First, the periodic potential which certainly induces a rapid 
density variation puts certain boundary conditions on the non-interacting 
energy eigenfunctions: that they be Bloch functions. As discussed in the 

Introduction, there still is a conserved momentum- like quantity associated 
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with this boundary condition, the crystal momentum. Thus it is possible 
to Fourier transform this many- body problem to momentum space pro- 
vided one remembers to allow for inter-band as well as intra-band 
transitions in any scattering events. Crystal momentum conserving 
delta functions appear and the problem thus becomes computationally 
manageable . 

Similar things happen in the impurity problem. The existence of an 
impurity potential which is singular at the origin requires all the single 
electron wavefunctions to readjust themselves, such that they have zero 
amplitude at the impurity. The self-consistent, many-body formalism of 
this report was just a fancy way of seeing to it that this boundary condition 
was satisfied. 

In the case of the model surface problem, the infinite barrier places 
the zero amplitude boundary condition on non-interacting electron wave- 
functions. In the cases considered here, sine waves satisfy this boundary 
condition. But sine waves are zero momentum eigenfunctions, so once 
again we see that there is something we can say about momentum con- 
servation in this non- invariant system. As has been seen, this fact 
enables us to go to momentum space where we can pick up a delta function, 
thus making calculations feasible. 

The following observations suggest a general rule. In all cases of 
rapid density variation, the cause has been some natural condition which 
imposes boundary conditions on s ingle- electron wavefunctions . The density 
variation, and thus long-range Friedel density oscillation, whether around 
an impurity or surface, is just a manifestation of the non-interacting 
electrons interfering in a manner which washes out the effects of the 
boundary conditions in a distance of order of Fermi wavelengths. The 
Friedel oscillations are not caused by any strange many-body effect, but 
only by a classical wave diffraction effect. In any event, for most rapidly 
varying density phenomena in nature, the cause can be seen to be a 
boundary condition which all electrons must satisfy, regardless of their 
energies. Perhaps this is the kernel of the idea: if one has a non- 

momentum or energy-dependent boundary condition, then one can usually 
obtain a set of single- electron basis functions satisfying this condition, 
for which something can be said about momentum eigenvalues. If this is 
the case, then it is possible to work in momentum space. As we have 
seen, for the lattice and surface problem, the eigenfunctions might be 
small sums of plane waves. Then the only problem when doing many-body 
theory in this basis is to keep track of interferences between various 
plane wave components of a single electron wave function. This might be 
a minor nuisance but some conservation laws are available for use. 

On the other hand, the slowly varying density variations considered 
by Kohn are usually man-made. The boundary conditions they impose are 
energy- dependent. Such things as applied fields, stresses, and minor 
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deformations on the atomic level are well treated within his formalism. 

In these sorts of problems it is realistic to think in terms of local 
microscopically uniform systems, along the philosophy of Thomas -Fermi 
theory, and develop a calculational scheme based on density- gradient 
expansions. But in these sorts of problems there is a one-to-one 
correspondence between the non-uniform, slowly varying potential pro- 
ducing the inhomogeneous gas and the non-uniform density, whereas in 
the diffraction- produced inhomogeneous electron gases, such is not the 
case. Thus it seems that the slowly and rapidly varying inhomogeneous 
electron gases are of fundamentally different origins and thus should be 
considered from fundamentally different theoretical viewpoints. 

The rule of thumb is that in a rapidly varying system, there is a 
significant boundary condition to be faced. Somehow one should get a 
reasonable single- electron basis set satisfying this boundary condition. 
There will always be something definite and simple that can be said 
about the momentum spectrum of this non-interacting basis. Then one 
should go to momentum space and keep track of all many- body interaction 
interference terms since interference effects are at the heart of this type 
of inhomogeneous electron gas. For the slowly varying case, one should 
forget momentum conservation and work within the density- gradient ex- 
pansion formalism of Kohn. Fortunately, most of the physically interesting 
systems fall into one or the other of these categorizations. 

The author feels that the realization of this point of view and the 
illustration of it by several specific examples constitute the major achieve- 
ment and positive contribution of this report. The specific conclusions of 
each illustrative example have been stated in their respective sections so 
there is little need to repeat them here. 

From a point of apparent chaos, the problem of many- body effects in 
rapidly varying inhomogeneous electron gases has been reduced to one 
which is conceptually little more difficult than the case of homogeneous 
electron gases, although a bit more troublesome from the computational 
point of view. This is because interferences between various , components 
of a single electron wavefunction must be properly catalogued. 

It is the hope of the author that the procedures developed and expounded 
here will be of use in future refinements of electron gas theory which will 
ultimately lead to the accurate, quantitative predictions of observable 
properties of real metals. 
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